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ABSTR4CT 


A  wind  tonnel  imestigatloa  vas  nade  at  Haeh 
muidbers  3  and  ^  of  the  heat  transfer  and  pressure 
distribution  on  a  flat  plate  in  a  region  where  an 
oblique  shockwave  Inqiinges  upon  and  Interacts  with 
a  turbulent  boundary  layer*  The  strength  of  the 
shockwave  was  varied  so  that  Interaction  regions 
both  with  and  without  separation  were  studied. 

Plots  of  the  ratio  of  the  local  Stanton  number 
through  the  interaction  region  to  the  zero 
pressure  gradient  Stanton  number  gave  qualitatively 
similar  distributions  for  both  Mach  numbers  and  all 
shock  strengths.  The  data  from  the  present  tests 
plus  the  results  of  e:q>eriments  reported  in  the 
literature  were  enq>loyed  in  the  construction  of 
a  method  which  predicts  the  heat  transfer 
distribution  In  an  Incident  shockwave  turbulent 
boundaiy  layer  interaction  region  over  a  fAlrly 
wide  range  of  free-stream  conditions* 
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sumuor 


A  wind  tunnel  Innreetigatlon  vae  made  of  the  leothermal  vail  heat 
transfer  rates  and  static  pressure  distributions  associated  vlth 
the  Intersctlffli  of  an  oblique  shockwave  vlth  a  turbulent  boundarjr 
lajer  on  a  flat  plate.  The  tests  were  conducted  at  nominal  Midi 
nuiDbers  of  3  and  $  and  free-strean  ReTnolds  nuniDersper  foot  of 
7«U  X  10^  and  3.6  x  ICP^  respeetlveljr.  All  measurements  were  made 
on  a  flat  test  plate  which  contained  Instrumentation  for  determining 
heat  transfer  rates,  wall  ten^ierature,  and  static  pressure  distribution.. 
Zn  addition,  two  boundary  layer  probes  were  permanently  Installed  on 
the  test  plate  at  separate  stations.  At  Kaeh  3  natural  transition 
of  the  test  plate  bound  azy  layer  occurred  upstream  of  the  area 
instrumented  for  heat  transfer  measurements.  A  system  of  injecting 
air  into  the  boundazy  layer  near  the  leading  edge  of  the  test  plate 
was  en^loyed  to  thicken  the  boundary  layer  at  Kach  3.  At  Hadi  $ 
it  was  foimd  necessary  to  e2i^}loy  both  grit  and  air  injection  into 
'  the  boundazy  layer  in  order  to  obtain  turbulent  flow  upstream  of  the 
heat.trazisfer  l^tzumentatlon.  Blowizig  an  additional  amount  of  air 
into  the  boundary  layer  above  the  minimum  amount  required  to  insure 
turbulence  did  not  appreciably  thicken  the  boundazy  layer  at  Kach  5. 

The  oblique  shockwave  which  inpinged  upon  the  test  plate  was  generated 
by  a  deflected  flat  plate  placed  above  the  test  plate.  The  angle 
of  the  shock  generator  plate  with  respect  to  the  wind  tunnel  center¬ 
line,  fiC  »  could  be  varied  from  0^  to  At  both  Mach  3  and  $ 
tuzuiel  uhokihg  occurred  at  Ad  somewhat  greater  than  12.0”. 

Both  the  test  plate  and  the  shock  generator  plate  spanned  the  wind 
tunzisl  test  section.  Observations  of  the  streamline  pattern  on  the 
test  plate  with  flow  paint  indicated  that  the  flow  field  was  two- 
dimensional  in  the  region  in  which  heat  transfer  and  pressure  measure- 
Bsnts  wexe  made. 

The  length,  ,  of  the  oblique  shockwave  turbulent  botoidary 
interaotlon  region  was  determined  from  the  static  pressure  distribu¬ 
tion.  It  was  found  that  the  rise  in  the  isothermal  heat  transfer 
coefficient,  hr  $  upstream  undisturbed  value,  , 

to  a  peak  value,  ,  occurred  well  within  the  pressure  interaction 
region  azid  covered  a  length  of  approximately  0.60  .  The  decay 

in  heat  transfer  coefficient  downstream  of  the  peak  value  was  fairly 
^slow  when  the  aft  end  of  the  test  plate  was  not  unduly  disturbed 
by  the  eaq;>ansion  from  the  trailing  edge  of  the  shock  generator  plate. 
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Attenpta  at  oorrelating  the  data  vara  nada  vslng  tha  local  Stanton 
miribarf  ,  and  the  tranafornad  monsatoa  thleloieaa  6^  as 
dafined  tgr  Raahotko  and  Taoker  In  NACA  IN  Ul5U* 

Plots  of  tha  Stanton  unnbor  ratio,  «  (anbaeript  d.  rafara 

to  local  eondltiona  and  /  to  tqpstreaa  nndlaturbad  conditions) 
la  tama  of  a  non-dtmenalonal  leni^  paranstar  revaaled  that  an 
apprcbctmata  alntlarlty  exlatad  among  all  of  tha  ShoelonTa  torbaleot 
bowdary  layer  data  obtained  in  tha  present  study.  Results  of  other 
published  experinenta  eonoamed  with  oonsiderably  dlffarent  nodals, 
bat  still  involving  the  intaraotion  of  a  turbulent  boundary  layer 
with  a  shockHava,  revealed  dlatribntions  of  vhleh  ware 

qualitatively  siMlar  to  the  present  tests.  All  the  suitable  available 
data  vara  used  to  derive  a  tentative  relationship  between  the  peak 
value  of  the  Stanton  nanber  ratio  and  the  static  pressure  ratio 
across  an  interaotion  region.  This  relationship  was  enployed  in 
the  eonstruotion  of  a  ealeulation  design  procedure  which  covers 
a  fairly  vide  range  of  free-strean  ecadltlons. 

Another  correlation  appeared  feasible  which  involved  the  transfomed 
Bonentum  thickness  raUo  across  the  heat  transfer  interaction  region. 
The  values  of  the  ratio  fall  Into  rou^ly  three  groups  characterized 
by  attached,  incipient,  and  separated  flow  conditions.  The 
correlation  indicated  that  the  R^olds  analogy  factor  decreases 
across  the  interaotion  region.  This  is  in  qualitative  agreement 
vlth  the  observed  decrease  in  the  local  recovery  factor  across  the 
interaotion. 

A  oorrelatlon  of  the  pressure  interaction  length  indicated  that 
whenever  air  Injection  thickened  the  boundary  layer,  the  results 
ware  equivalent  to  a  boundary  layer  thickened  by  an  extension  of 
tbs  boundary  layer  run.  The  heat  transfer  measurements  were 
practically  Insensitive  to  air  Injectim. 
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I.  INTRQDOCTION 

On  a  vehicle  traveling  at  supersonic  speeds  there  will  usually  be 
found  one  or  more  shockwaves  which  impinge  upon  or  emanate  from 
a  region  in  which  a  boundary  layer  has  developed.  In  cases  of 
supersonic  air-breathing  propulsion  systems  the  occurrence  of 
slwckwave  Impingement  on  the  Inlet  diffuser  boundary  layer  is 
unavoidable*  The  interaction  of  a  shockwave  with  a  boundary  layer 
will  modify  both  the.  viscid  and  invlscid  flow  field.  Inviscld 
shockwave  theory  specifies  a  step  function  pressure  distribution 
across  the  shock  front.  However,  experimental  evidence  shows  that 
the  rise  in  pressure  across  the  shockwave  is  obtained  over  a  finite 
length  in  the  region  of  interaction  with  the  boundary  layer.  Of 
greater  significance  are  the  results  of  a  few  tests  which  indicate 
that  the  heat  transfer  rates  in  the  aft  portion  of  the  interaction 
are  higher  than  would  be  computed  from  undisturbed  local  flow 
conditions.  Since  the  temperatures  encountered  at  flight  speeds 
in  excess  of  Mach  3  may  seriously  affect  a  vehicle’s  structural 
integrity  it  is  clear  that  an  accurate  prediction  of  heat  transfer 
rates  in  an  area  influenced  by  shockv;ave  boundary  layer  interaction 
is  important. 

Of  course,  the  interaction  of  a  shockwave  with  a  boundary  layer 
is  a  highly  complex  phenomenon.  No  theory  has  been  devised,  even 
for  restricted  conditions,  which  correctly  predicts  the  effect 
of  the  interaction  on  all  of  the  fluid  mechanical  properties. 
Therefore,  the  necessity  of  obtaining  experimental  data  in  this 
field  is  obvious.  It  was  found  by  Chapman  (Reference  9)  and  others 
that  it  is  meaningful  to  divide  all  shockwave  boundary  interactions 
into  three  basic  types,  viz.,  (1)  the  boundary  layer  is  laminar 
throughout  the  interaction,  (2)  the  boundary  layer  is  laminar  at  the 
beginning  of  the  interaction  but  transition  occurs  before  the  end 
of  the  interaction,  and  (3)  the  boundary  layer  is  turbulent  upstream 
of  the  interaction  and  remains  turbulent  throughout  the  interaction. 
The  present  study  is  primarily  concerned  with  the  heat  transfer  rates 
associated  with  the  third  type  of  interaction.  Unfortunately, 
almost  all  of  the  shockwave  turbulent  boundary  layer  interaction 
experiments  found  in  the  literature  in  which  a  useful  range  of 
parameters  were  tested  contain  only  pressure  measurement  data 
(e.g.  see  References  1,  2  and  9).  The  limited  amount  of  available 
experimental  heat  transfer  data  which  pertains  to  the  present 
problem  falls  into  two  classes.  First,  the  tests  were  made  on 
specific  vehicle  configurations,  with  the  resulting  flow  field  so 
complex,  that  analysis  or  correlations  with  basic  variables  is 
virtually  Impossible.  In  the  second  class  of  experiments,  tests 
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have  been  made  in  a  twoodinensional  flow  field,  bat  with  only  one 
Maeh  number  and  one  or  two  ehock  strengths  so  that  although  the 
flow  field  can  be  defined  the  data  range  is  insufficient  for 
correlation.  In  this  latter  category  are  the  experiments  described 
in  References  (12),  (15)>  (16)  and  (17).  In  these  tests  the 
shockware  Is  generated  by  the  surface  along  which  the  boundary  layer 
has  developed .  No  case  was  found  in  which  the  heat  transfer  data 
were  obtained  in  the  region  of  an  oblique  shockwave  impinging  upon 
a  turbulent  boundary  layer. 

The  present  tests  were  conducted  in  order  to  partially  fill  the 
considerably  large  gap  in  knowledge  concerning  the  experimental 
heat  transfer  associated  with  shockwave  turbulent  boundary  layer 
interaction.  In  particular,  a  two-dimensional  experimental  heat 
transfer  Investigation  was  made  at  Mach  numbers  3  and  5  on  a  flat 
plate  with  a  turbulent  boundary  layer  and  with  an  impinging  oblique 
shockwave  of  variable  strength.  An  attempt  was  made  to  increase 
the  boundary  layer  thickness  utilizing  an  air  injection  system. 

A  sufficient  amount  of  data  was  sought  so  that  the  heat  transfer 
rates  could  be  predicted  in  the  interaction  region  over  a  useful 
Mach  number  and  Reynolds  number  range. 
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II.  EXPERIMENTAL  PROGRAM 


A.  EXPERDffiNTAL  APPARATUS 

The  model  vas  designed  to  provide  the  following: 

(a)  Two-dimensional  flow  on  a  plane  surface 
(the  heat  transfer  test  plate), 

(b)  An  oblique  shock  of  variable  strength  impinging 
on  the  boundary  layer  on  the  test  plate. 

(c)  A  boundary  layer  on  the  test  plate  which  becomes 
turbulent  upstream  of  the  influence  of  the 
impinging  shock, 

(d)  Instruniontation  on  the  test  plate  for  determining 
both  the  heat  transfer  coefficient  and  the  static 
pressure  distribution  in  the  shockwave  boundary 
layer  interaction  region  and  the  boundary  layer 
velocity  profiles  upstream  and  downstream  of  the 
interaction. 

Schematics  of  the  model  and  certain  details  of  the  Instrumentation 
are  shown  in  Figures  1  through  7.  Photographs  of  the  model  are  shown 
in  Figures  8  and  9* 

The  results  of  the  experiments  described  in  References  (l)  and  (2) 
indicate  that  two-dimensional  flow  and  shock  systems  can  be  obtained 
only  if  the  tost  plate  and  the  shock  generator  plate  are  essentially 
free  of  aspect  ratio  effects.  In  the  present  caso  this  was 
accomplished  by  designing  the  model  so  that  it  completely  spanned 
the  tunnel  test  cccticn.  Both  the  test  plate  and  the  shock  generator 
plate  were  supported  from  circular  steel  plates  which  replaced  the 
existing  Schliercii  windows  in  the  tunnel  (Figure  1).  The  test 
plate  attitude  could  be  adjusted  within  the  limits  of  +2°  in  order 
to  permit  alignment  with  the  free  stream.  The  shock  generator  plate 
attitude  could  be  varied  from  0®  to  15®.  The  generator  plate  had 
a  different  pivot  point  for  ^ach  frcr-strc~r  Mr.ch  number  which 
allowed  the  shock  to  impinge  approximately  11.75"  from  the  leading 
edge  of  the  test  plate  (i.e.,  the  window  centerline)  for  generator 
plate  angles  from  0  to  12®.  Angles  above  12®  caused  tunnel  choking. 
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The  test  plate  was  provided  with  a  botmdary  layer  air  Injection 
eyetem  which  was  used  for  either  thickening  an  already  turbulent 
boundary  layer  or  helping  to  Induce  full  turbulence  In  a  transitional 
boundary  layer.  The  air  Injection  system  (see  Figure  2)  consisted 
of  a  row  of  holes  2.58  Inches  downstream  of  the  leading  edge  of  the 
test  plate,  a  plenum  chamber  tindemeath  the  plate,  and  the  necessary 
plumbing  and  control  valves  to  provide  a  controllable  amount  of 
airflow  Injected  normal  to  the  free-stream  airflow.  FOr  the  Mach  5*0 
run  a  wide  fixed  transition  strip  of  .033”  grit  was  located  1.5” 
(model  station  1.5)  from  the  leading  edge  of  the  test  plate. 

The  static  pressure  distribution  on  the  test  plate  was  measured  with 
a  row  of  twenty-three,  .05l”  Inside  diameter,  flush,  pressure  taps 
located  2.5”  off  the  model  centerline  and  placed  apart  with  the 
first  tap  located  at  model  station  3.5  (see  Figure  U).  Pitot  rakes 
were  permanently  Installed  to  measure  the  boundary  layer  velocity 
profiles  at  modal  stations  9.0  and  15.0.  Pressure  measurements 
were  made  In  the  boundary  layer  air  injection  system  to  pennlt 
duplication  of  blowing  rates.  One  pressure  was  measured  In  the 
air  supply  tube  between  a  needle  valve  and  a  restrictor  located  outside 
of  the  tunnel.  A  second  pressure,  ,  was  measured  In  the  plenum 
chamber  \inder  the  leading  edge  of  the  test  plate. 

The  heat  transfer  Instrumentation  on  the  test  plate  Is  similar  to 
that  used  In  References  (3)  and  (U).  A  single  heat  transfer  gage 
consists  of  a  nlchrome  ribbon  heating  element  cemented  to  an 
Insulating  surface  (melamine  fiber  glass  composite)  plus  a 
thermocouple  which  measures  the  element  temperature.  The  thermocouple 
was  inserted  In  the  fiber  glass  approximately  .003"  below  the  center 
of  the  element.  Since  fiber  glass  is  an  electrical  as  well  as  a 
thermal  Insulating  plastic,  the  thermocouple  was  electrically  Isolated 
from  the  nlchrome  ribbon.  The  electrical  current  flow  path  in  the 
heating  elements  was  made  by  milling  slots  (.002"  to  .003"  wide) 

In  the  metal  as  shown  In  Figure  2.  Originally  the  slots  were  filled 
with  Krylon  in  order  to  obtain  a  smooth  surface.  During  the  Initial 
phase  of  the  wind  tunnel  tests,  however.  It  was  found  that  the  Krylon 
was  bubbling  and  peeling  off.  Attempts  to  use  other  Insulating 
filler  material  had  the  same  undesirable  results.  After  a  discussion 
with  an  experimenter  who  had  tested  a  similar  heat  transfer  model. 

It  was  concluded  that  the  thickness  of  1he  nlchrome  ribbon  (about 
.002")  would  have  a  negligible  effect  on  the  results  since  this 
dimension  represented  less  than  2%  of  the  boundary  layer  thickness. 
Consequently,  the  slots  were  left  unfilled  for  the  remaining  runs. 

Figure  5  Is  a  schematic  of  a  typical  heater  element  circuit.  It 
is  seen  that  a  variable  resistor  was  used  to  control  the  current 
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in  the  heater  element.  The  current  was  supplied  hy  a  constant 
Toltage,  direct  current  power  supply. 

Spanwlse  heat  loss  from  the  transverse  heater  elements  was  minimized 
by  guard  heater  elements  located  along  the  sides  of  the  heat  transfer 
area  (see  Figure  li).  Four  thermocouples  were  located  on  the  back  of 
'Uie  fiber  glass  slab  to  aid  in  determining  the  conductive  heat  loss 
through  the  fiber  glass  slab.  Also,  a  thermocouple  was  inserted  in 
the  lower  surface  of  the  shock  generator  plate  to  measure  the  recovery 
tenqperature  for  use  in  calculating  the  radiation  heat  loss.  In  all, 

38  heating  elements  including  the  guard  elements  and  3^  thermocouples 
were  installed  on  the  test  plate,  with  an  additional  thermocouple 
on  the  shock  generator. 

At  least  one  Schlleren  photograph  was  taken  for  each  run  en^loying 
the  windows  located  In  the  steel  plates  which  supp>orted  the  model. 
These  windows  were  only  three  inches  in  diameter  and  thus,  only  a 
small  portion  of  the  flow  field  could  be  viewed.  Larger  windows 
were  not  used  because  cf  various  design  considerations.  The 
Schlleren  system  at  the  test  facility  was  designed  to  employ  a 
collimated  light  beam  approximately  I8  inches  in  diameter.  Since 
the  system  was  basically  Incompatible  with  the  three  inch  window 
diameter,  it  is  not  too  surprising  that  the  photographs  suffer  from 
poor  definition.  Schlleren  photographs  at  Mach  3.0  corresponding 
to  each  angle  of  the  shock  generator  plate  tested  are  shown  in 
Figures  10  through  lit.  Only  one  photograph.  Figure  15,  taken  at 
Hach  is  shawn.  At  Mach  5.0  almost  all  the  interaction  occurred 
upstream  of  the  region  viewed  by  the  photographs. 


B.  EXPERIMENTAL  PROGRAM  AND  DATA  REDUCTION 

The  experimental  program  was  conducted  at  the  USN  Missile  and 
Astronautics  Test  Center,  Point  Mugu,  California,  in  the  time 
period  from  7  May  through  25  May  1962.  The  tests  were  conducted 
at  free-stream  Mach  numbers  of  2.95  and  5.02  (nominally  3.0  and  5.0) 
and  at  Itee-streara  Reynolds  numbers  per  foot  of  7.U  x  10°  and 
3.6  X  10®,  respectively.  A  summary  of  the  test  conditions  at 
each  Mhch  number  is  given  in  Table  I. 

The  heat  transfer  measurements  on  the  test  plate  were  limited  to 
2h  data  channels  by  the  equipment  available  at  the  test  facility. 
That  is,  only  the  tenperature  and  electrical  power  dissipation  of 
2h  heater  elements  could  be  monitored  for  a  given  heat  transfer 
test.  The  side  guard  heater  elements  had  to  be  included  in  the 
2U  heater  element  limitation.  Therefore,  several  preliminary 
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test  runs  were  made  at  each  Mach  number  in  order  to  establlah  the 
shoelwave  boundary  layer  Interaction  length  which  in  turn  set  a 
lower  limit  on  the  allowable  length  of  the  instrumented  section 
to  be  used.  The  heater  element  patterns  selected  for  Kach  numbers 
3  and  $  are  shown  in  figures  6  and  7>  respectively.  It  may  be 
noted  by  comparing  Figures  6  and  7  with  Figure  U  that  some  of  the 
vide  nichrome  ribbons  are  connected  in  series  to  operate  effectively 
as  a  I"  wide  element.  Although  the  output  of  all  thirty-six 
thermocouples  were  read  into  the  data  system  for  each  data  point, 
only  selected  thermocouples  were  used  for  heat  transfer  data 
reduction. 

A  heat  transfer  test  run  consisted  of  data  points  taken  at  three 
or  fo\ir  different  temperatures  with  the  shock  strength  and  the 
condition  of  the  boundary  layer  upstream  of  the  Interaction  held 
constant.  A  data  point  was  obtained  by  adjusting  the  current  flow 
In  each  heater  element  (including  the  guard  elements)  so  that  all 
thermocouples  used  for  monitoring  the  elements  read  essentially  the 
same  temperature  .  The  heat  transfer  coefficient,  ^  ,  measured  by 
a  given  element  was  obtained  by  plotting  the  convective  heat  loss 
per  unit  area  per  second,  C^/a)c  »  versus  the  temperature  of  the 
nichrome  ribbon,  77  (see  Figure  16).  The  slope  of  the  straight  line 
(determined  by  a  least  square  fit)  through  these  points  yields  the 
heat  transfer  coefficient  through  ‘Uie  relationship. 


where,  as  shown  in  Figure  16,  the  effective  recovery  temperature,  %  , 
is  defined  as  the  intercept  of  the  straight  line  on  the 
versus  %  plot  with  the  line  o  .  The  values  of 

and  are  the  result  of  corrections  to  the  actual  measurements. 

The  corrected  convective  heat  loss  is  given  as 


(xl- 


’  power  conducted 
and  radiated 

(electrical  power  dissipated  by  the  element )-\ away  from  the  element 
pianform  area  of  the  element 


The  heat  loss  by  conduction  was  calculated  by  assuming  a  one -dimensional 
heat  flow  to  the  back  of  the  fiber  glass.  The  thermocouples  beneath  the 
guard  elements  indicated  that  this  as5iu)g>tion  was  accurate.  However, 
when  a  shock  impinged  on  the  test  plate,  at  least  one  pair  of  side 
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guard  elements  was  in  a  region  of  strongly  varying  heat  transfer. 

With  the  aid  of  a  digital  computer  program  calculations  were  made 
of  the  spanwise  heat  loss  from  the  heat  transfer  elements  through 
the  guard  elements  in  a  region  of  varying  heat  transfer.  It  was 
found  that  for  the  worst  case  at  Mach  $  the  conduction  losses  through 
the  guard  elements  represented  less  than  h%  of  the  heat  input  to  a 
heat  transfer  element.  At  Mach  3  the  effect  would  be  even  less 
significant.  The  radiation  correction  was  obtained  with  the  aid 
of  the  measured  temperatures  of  the  wind  tunnel  wall  and  shock 
generator  plate  and  employed  the  following  assumptions: 

(a)  The  emissivity  of  the  nichroma  ribbon  was  0.7. 

(b)  The  view  factor  from  a  heater  element  to  the  shock 
generator  plate  was  0.87. 

(c)  The  view  factor  from  a  heater  element  to  the  wind 
tunnel  walls  was  0.13. 

For  the  undisturbed  flat  plate  the  combination  of  conduction  and 
radiation  losses  represented  approximately  10^  of  the  electrical 
energy  dissipated  in  an  element  at  Mach  3.0  and  UO^  at  Mach  5.0. 

Although  the  thermocouple  employed  for  monitoring  the  temperature 
of  a  particular  heater  element  was  located  only  .003"  beneath  the 
element,  it  is  not  correct  to  equate  the  temperature  of  the 
element,  7i/  »  with  the  thermocouple  temperature,  T  .  This  is 
because  of  the  large  difference  in  thermal  conductivity  of  the 
thermocouple  as  compared  with  the  surrounding  fiber  glass  and  the 
resulting  distortion  of  the  heat  flow  field.  An  investigation  was 
made  employing  a  digital  computer  program  to  determine  the  value 
of  7}if  for  a  given  value  of  T  .  It  was  found  that  7?  could 
be  expressed  approximately  in  the  form  where  «.  and  A 

are  constants  and  T,  is  the  temperature  recorded  on  the  back  of 
the  fiber  glass  slab.  The  highest  value  of  T  used  for  most  of 
the  runs  was  about  160®F  which  yields  a  value  of  165®F  for  . 

This  test  was  conducted  so  that  the  heat  transfer  coefficients 
could  be  determined  for  the  case  of  an  isothermal  wall.  Thus, 
the  effect  of  deviations  from  this  ideal  condition  must  be  evaluated 
in  order  that  proper  corrections  can  be  made  to  the  experimental 
results.  Within  the  area  instrumented  for  heat  transfer  measurements 
and  surrounded  by  heater  guard  elements,  the  variation  in  temperature 
was  small  enough  so  that,  in  general,  no  correction  was  necessary. 
However,  since  the  heated  portion  of  the  test  plate  did  not  extend 
to  the  leading  edge,  an  important  correction  had  to  be  made  due  to 
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the  step  function  rise  In  temperature  and  its  consequential  effect 
on  the  experimentally  determined  heat  transfer  rates.  In  Reference  (5) 
the  effect  of  a  step- tempera tore  distribution  on  the  heat  transfer 
coefficient  for  turbulent  flow  over  a  flat  plate  with  zero  pressure 
gradient  is  given  as 


where 


/  - 


heat  transfer  coefficient  for  step  temperature 
distribution 


heat  transfer  coefficient  for  a  true  isothermal  wall 


V  ■  reciprocal  of  exponent  in  power  law  for  boundary 
layer  velocity  profile 

L  ■  distance  from  leading  edge  of  the  plate  to  the 
ten^rature  step 


'K  ■  distance  from  leading  edge 

K  Is  a  function  of  ^4^  ,  Mach  number,  and  }  however, 

for  values  of  near  \mity  (as  in  this  test)  a  value  of  AT’-/ 

is  correct.  For  the  present  experiments  can  be  considered 

as  a  typical  value  valid  for  both  the  Mach  3  and  $  tests.  For  the 
case  of  shockwave  boundary  layer  interaction  the  correction  was 
assumed  to  be  the  same  as  for  a  flat  plate  with  zero  pressure 
gradient.  Therefore,  all  the  heat  transfer  data  were  corrected 
to  Isothermal  wall  conditions  by  the  equation. 


i 


where  kc  is  the  heat  transfer  coefficient  with  the  corrections  made 
for  conduction,  radiation,  and  thermocouple  reading. 
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It  Is  to  be  noted  that  at  each  Mach  number  almost  all  the  runs  were 
made  at  one  of  two  values  of  the  boundary  layer  air  Injection  plenum 
chamber  pressure  ratio,  .  At  Mach  3  the  two  values  were 

■  1.0  and  5.9.  •  1*0,  which  corresponds  to  no  boundary 

layer  blowing,  was  normally  used  at  Mach  3  since  the  boundary  layer 
was  turbulent  upstream  of  the  heat  transfer  instrumentation.  For 

>  6  it  was  found  that  the  boundary  layer  velocity  profiles  became 
distoirbed  and  thus  blowing  was  producing  an  undesirable  effect.  Since 
the  runs  at  ■  5.9  did  not  produce  large  changes  in  the  static 

pressure  and  heat  transfer  distribution  through  an  interaction  region, 
just  two  exploratory  runs  at  an  intermediate  value  ■  2.0) 

were  made. 

At  Mach  5.0  it  was  necessary  to  employ  boundary  layer  blowing,  in 
addition  to  the  wide  strip  of  grit,  merely  to  insure  turbulence 
in  the  boundary  layer.  The  minimum  value  of  which  gave 

turbulent  flow  was  approximately  3.5.  It  is  to  be  noted  that  ^/f>^  ^  l.U 

corresponds  to  no  boundary  layer  blowing  at  Mach  5  because  the  static 
pressure  on  the  undisturbed  test  plate  was  about  l.U  times  the 
tabulated  free-stream  pressure.  For  "y  10.0  the  velocity  profiles 

became  distorted.  In  general,  blowing  had  a  smaller  effect  at 
Mach  5  than  at  Mach  3*  therefore,  almost  all  the  turbulent  data 
were  obtained  at  blowing  rates  iii  the  vicinity  of  the  two  allowable 
extremes  of  . 

Several  difficulties  arose  in  interpreting  the  Mach  5  tabulated 
data.  During  the  Mach  5  tests  it  was  noticed  that  thermocouple  #18 
(see  Figure  U)  beneath  the  heater  element  #13  (see  Figure  7) 
recoixled  a  temperature  considerably  lower  than  the  thermocouples 
placed  directly  below  the  other  nichrome  elements.  At  the  time, 
this  anomaly  was  thought  to  be  due  to  a  malfunction  in  the  thermocouple. 
It  was  not  considered  serious  because  another  thermocouple  (#17) 
was  used  to  monitor  the  temperature  of  heater  element  #13.  After 
the  tests  were  completed,  it  was  found  that  the  thermocouple  #18 
had  not  malfunctioned.  Upon  disassembling  the  test  plate  instrumenta¬ 
tion  it  was  found  that  the  wide  nichrome  strip-placed  directly 
above  thermocouple  #18  (this  nichrome  strip  was  the  downstream  half 
of  the  element  A3)  was  subjected  to  a  partial  electrical  short. 

That  is,  a  portion  of  the  recorded  electrical  current  found  a  path 
through  the  steel  plate  rather  than  through  the  nichrome  strip. 
Calculations  were  made  utilizing  the  recorded  temperatures  in  the 
vicinity  of  the  partial  electrical  short  in  order  to  obtain  an 
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upper  bound  on  Its  effect  on  the  calculated  heat  transfer  coefficients. 
It  vas  found  that  a  sizeable  error  in  ^  due  to  the  electrical 
short  could  have  occurred  only  for  heater  elements  #13  and  #IU.  This 
error  is  no  larger  than  30%  (and  probably  much  less)  for  a  zero  pressure 
gradient.  Vithod?  elements  #13  and  #lli  are  in  the  interaction 
region,  therefore,  the  error  is  not  as  significant,  being  less  than 
B%  for  fiC  ■  Since  the  actual  percentage  coirrection  cannot  be 

determined  and  since  conclusions  drawn  from  the  data  are  not  seriously 
affected,  /r  was  not  corrected  for  the  error  arising  from  the 
partial  electrical  short. 

A  further  difficulty  arose  with  respect  to  interpreting  the  Mach  $ 
data  and  by  coincidence  heater  elements  #13  and  #lU  were  .involved. 

The  tabulated  data  from  Pt.  Mugu  indicated  that  for  <.>  0°  the 
value  of  /ir  measured  with  heater  elements  #13  and  #lU  (the  center 
cf  each  element  is  at  model  stations  12.75  and  13.75,  respectively) 
vas  consistently  about  $0%  higher  than  /jr  measured  with  heater 
elements  #12,  #15,  and  #16.  This  "blip"  in  the  distribution  of  the 
tabulated  hr  could  not  be  attributed  to  the  expansion  wave 
emanating  from  the  trailing  edge  of  the  shock  generator  plate  since 
the  "blip"  occurred  even  when  the  pressure  distribution  was  relatively 
flat  downstream  of  heater  element  #12.  It  is  felt  that  the  "blip"  is 
not  the  result  of  anomalous  behavior  of  the  flow,  but  rather  it  is 
due  to  an  error  in  the  data  reduction  procedure.  Specifically,  the 
error  consisted  of  assigning  the  measured  electrical  current  in  the 
guard  elements  #23  and  #2U  to  #13  and  #14  and  vice  versa.  Whether  the 
nisarrangement  was  in  the  electrical  connections  to  the  data  reduction 
system  or  was  within  the  data  reduction  program  itself  is  not  known. 
Meither  possibility  can  be  checked  since  the  model  has  been  dismantled 
and  the  Pt.  Mugu  facilities  have  been  closed.  Nevertheless,  self- 
consistency  of  the  Mach  5  data  strongly  implied  that  the  error 
actually  occurred  and  the  data  were  corrected  with  this  in  mind. 

Figures  47  through  60  which  give  the  Mach  5  heat  transfer  coefficient 
distributions  also  contain  the  data  points  (flagged  symbols)  that 
are  felt  to  be  erroneous.  These  "incorrect"  data  are  shown  because, 

(l)  they  are  contained  in  the  original  data  reduction  and  tabulation 
i^e  at  Pt.  Mugu  and  (2)  It  is  not  possible  to  assign  a  100^ 
certainty  to  the  belief  that  the  flagged  data  points  are  truly 
incorrect. 
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III.  EXPERIMENTAL  RESULTS 


In  this  section  the  Kach  3.0  and  Mach  $.0  data  are  presented 
separately.  At  each  Mach  number  the  most  significant  results 
are  the  distributions  of  static  pressure  on  the  test  plate,  corrected 
heat  transfer  coefficient,  and  recovery  temperature.  Also,  the 
velocity  profiles  obtained  at  model  stations  9.0  and  1$.0  are  shown. 
Although  this  report  Is  only  concerned  with  turbulent  Interactions 
(i.e.,  the  boundary  la^r  Is  turbulent  upstream  of  the  Influence 
of  the  Impinging  shocks  supplemental  data  are  presented  which  were 
acquired  during  tests  In  which  the  boundary  layer  became  turbulent 
in  the  interaction  region. 


A.  MACH  3.0  DATA 

A  free-stream  Mach  number  of  2.9?  was  calculated  from  the  tunnel 
settling  chamber  pressure,  ^  ,  and  the  pressure  measured  with  a 
pitot  tube  placed  below  and  forward  of  the  test  plate.  The  static 
taps  on  the  test  plate,  however.  Indicated  pressures  somewhat  less 
than  the  free-stream  value  in  a  region  unaffected  by  the  interaction. 
If  the  test  plate  static  pressure  Is  assumed  to  be  the  result  of  an 
isentroplc  expansion  from  free-stream  conditions,  then  the  Mach 
number  in  the  undisturbed  region  over  the  test  plate  Is  2.98.  Thus, 
2.98  Is  considered  to  be  the  free-stream  Kach  number  of  the  shock- 
wave  boundary  layer  Interaction. 


1.  Static  Pressures  on  the  Test  Plate 

The  static  pressure  distribution  on  the  test  plate  for  the  various 
shock  generator  plate  angles,  ,  and  for  the  various  boundary  layer 
blowing  ratios,  ,  aro  shown  in  Figures  17  through  19.  Noted 

on  the  graphs  are  three  significant  quantities,  namely, 

(a)  The  point  at  which  the  shock  would  strike  the 
boundary  layer  if  there  were  no  Interaction. 

(Determined  graoblcally  from  scale  drawings  of  the 
test  apparatus.) 


( 
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(b)  Ihs  point  at  idtleh  the  expansion  from  the  trailing 
edge  of  the  shock  generator  plate  Impinges  on  the 
boundaiy  layer  downstream  of  the  Interaction. 

(Determined  in  same  manner  as  (a)). 

(e)  The  Invlscld  pressure  rise  which  results  from  a 
reflected  shock  generated  by  a  wedge  at  angle 
and  upstream  Kach  number  of  2.98. 

all  cases  the  static  pres8\ire  rose  above  the  Inviscid  value 
and  then  appeared  to  decrease  before  the  expansion  from  the  aft 
end  of  the  generator  plate  influenced  the  pressure.  At  ^  «  12^ 
the  generator  plate  expansion  Impinged  at  ^  •  1^.3"  and,  thus, 

■ay  have  affected  the  last  two  pressure  taps.  Fbr  12®  the 
expansion  Impinged  downstream  of 16"  and  should  not  have 
influenced  the  pressure  measurements.  In  similar  experiments  where 
an  incident  oblque  shock  Interacted  with  a  boundary  layer.  It  was 
found  that'  the  pressure  in  the  reat-tachment  region  will  occasionally 
exceed  the  inviscid  value.  However,  in  the  present  case,  the 
consistency  and  the  size  of  the  pressure  "blip"  ijiplies  that  either 
a  small,  additional,  unexplained  compression  and  then  expansion  were 
generated  somewhere  in  the  flow  field  or  else  'the  reattaching  flow 
overcompressed  and  then  expanded  to  the  proper  downstream  pressure. 

It  may  be  noted  that  a  small  pre.'ssure  "blip"  occurred  for  “0° 
in  the  region  where  a  Mach  line  emanating  from  the  leading  edge  of 
the  shock  generator  plate  struck  the  test  plate  boundary  layer.  This 
"blip"  pjrobably  arose  from  either  the  pressure  field  associated 
with  the  slight  bluntness  of  the  leading  edge  of  'the  shock  generator 
pla'te  (diameter  ^  .01")  or  the  pressure  field  generated  above  the 
generator  plate  which  was  'then  fed  upstream  through  the  sidewall 
boundary  layer. 

The  generator  plate  angles  •  6®,  9.5®»  and  12*^  were  chosen  'to 
correspond  to  the  conditions  of  fully  attached  flow.  Incipient 
separation,  and  separa'ted  flow,  respectively.  Incipient  separation 
may  be  defined  to  occur,  for  Increasing  shock  strength,  at  the  value 
of  ^ich  corresponds  to  a  change  in  the  static  pressure 
distribution  from  one  with  one  Inflection  point  'to  a  distribution 
with  three  inflection  points.  Kuehn  in  Reference  (6)  developed  this 
criterion.  For  the  present  tests  the  pressure  'taps  were  spaced 
apai*t  and  did  not  yield  a  sufficient  definition  of  'the  pressure 
dis'tributlon  to  determine  the  occurrence  of  Incipient  separation. 
However,  the  Schlieren  photograph  for  oC  ■  9.5°  (see  Figure  13) 
did  indicate  a  small  region  of  separated  flow.  For  ■  12®  the 
pressure  distribution  showed  three  distinct  inflection  points  which 
definitely  indicated  separated  and  reat'tached  flow.  The  Schlieren 
photograph  in  Figure  lit  verifies  this  observation. 
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Air  Injection  (i«e.«  blowing)  into  the  boundary  layer  has  a 
oonparatively  small  effect  on  the  pressure  distribution.  As  shown 
in  Figure  19  the  length  of  the  interaction  region  does  increase 
slightly  with  blowing  and  the  peak  pressure  appears  to  decrease 
slightly.  Within  the  accuracy  of  the  instrumentation  blowing  has  no 
influence  on  the  static  pressure  upstream  of  the  interaction. 


2,  Boundary  Layer  Rakes 

Although  the  boundary  layer  rakes  at  model  stations  9.0  and  1^.0 
did  not  contain  a  sufficient  number  of  pitot  tubes  for  obtaining  an 
accurate  velocity  profile,  they  were  useful  in  indicating  certain 
gross  effects.  The  profiles  obtained  at  Mach  2.95  are  ^own  in 
Figures  20  and  21.  The  thickness  of  the  boundary  layer,  S  ,  was 
determined  from  log-log  plots  of  distance,  ^  ,  above  the  test 
plate  versus  velocity  ratio,  .  The  results  appear  to  be 

accurate  to  +  0.02  inches.  With  no  blowing  and  no  shock,  S  was 
approximately  0.13"  and  0.19"  at  model  stations  9.0  and  15.0, 
respectively.  It  was  found  that  for  >  6  the  velocity  profiles 

became  distorted  at  station  9.0  without  increasing  S  .  For 

"  2.0,  S  increased  to  0.16"  and  0.23"  at  the  upstream  and 

downstream  probes,  respectively.  With  a  blowing  ratio  of  5.9» 

S  •  0.19"  at  the  forward  probe. 

Values  of  the  momentum  thickness,  0  ,  were  calculated  from  the 
velocity  (or  Mach  number)  profiles  in  several  cases  by  assuming 
constant  total  temperature  throughout  the  boundary  layer  and 
Integrating  the  appropriate  function  by  plotting  the  integrand 
against  distance  above  the  plate  and  using  a  planimeter  to  find 
the  area  under  the  curve. 


3*  Heat  Transfer  Coefficient 

The  heat  transfer  coefficient  distributions  (corrected  to  an 
isothermal  wall)  for  Mach  3  are  shown  in  Figures  22  through  33. 

On  each  graph  three  theoretical  lines  are  drawn  which  were 
determined  for  zero  pressure  gradient  (i.e.,  no  shock)  turbulent 
flow.  These  lines  are  the  results  of  calculations  of  the  Reshotko 
and  Tucker  theory  (Reference  7)  assuming  the  boundary  layer 
turbulent  at  the  leading  edge  of  the  test  plate.  The  three  lines 
correspond  to  the  following: 
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(a)  Heat  transfer  distribution  calculated  strictly 
by  the  method  of  Reference  (7)* 

(b)  Heat  transfer  distribution  as  given  by  Reference  (7) 

but  incorporating  the  experimentally  determined 
momentum  thickness  at  the  upstream  boundazy  layer 
probe  for  •  1.0. 

(o)  Heat  transfer  distribution  as  given  by  Reference  (7) 
but  incorporating  the  experimentally  determined 
momentum  thickness  at  the  upstream  boundary  laysr 
probe  for  “  5.9. 


It  is  interesting  to  note  that  the  line  corresponding  to  the  strict 
use  of  the  theory  lies  between  the  two  that  employ  the  measured 
momentum  thickness,  O  .  The  measured  Q  's  at  station  9.0  are 
5.7  X  lO"'^  and  13.7  x  10”^  inches  for  ■  1.0  and  5.9,  respectively. 

The  change  is  about  a  factor  of  2.1z  in  &  but  this  infers  only  a  20^ 
variation  in  heat  transfer  according  to  Reshotko's  theory.  Thus, 
it  is  not  too  surprising  to  find  that  for  the  undisturbed  boundary 
layer  the  difference  in  the  experimental  heat  transfer  between 
blowing  and  no  blowing  is  essentially  hidden  by  the  experimental 
scatter. 

It  may  be  seen  that  a  bump  occurs  in  the  isothermal  heat  transfer 
coefficient,  >  distribution  between  stations  12.5  and  13.5  when 
the  shock  generator  plate  is  removed  and  when  the  generator  is 
installed  at  od  ■  0°  (see  Figures  22  through  25).  The  reason  for 
this  rise  in  is  not  understood.  When  the  shock  generator  is 
installed  at  <a6  ®  0°  there  is  a  small  pressure  disturbance  emanating 
from  the  leading  edge  of  the  plate  (this  was  discussed  in  Section  III  A  1), 
but  this  Increase  in  pressure  does  not  warrant  the  fairly  sizeable 
increase  in  hr  .  However,  with  toe  generator  plate  removed  no 
pressure  disturbance  is  found  and  toe  area  under  the  bump  is  reduced. 

Thus,  this  anomaly  in  the  Ar  distribution  may  be  due  to  causes 
not  easily  detected,  such  as  cross  flow  due  to  asymmetric  disturbance 
of  toe  sidewall  bou^ary  layer.  It  is  believed  that  when  strong 
disturbances  are  Introduced  in  the  flow  by  deflecting  the  generator 
plate,  this  phenomenon  will  not  Introduce  any  important  errors  in 
toe  results.  One  Indirect  verification  of  thp  argument  is  shown  in 
Figure  71  where  good  correlation  of  pressure  Interaction  length,  , 
is  found  with  toe  present  data  and  the  results  of  other  experiments. 
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The  length  of  the  pressure  interaction  region  ascertained  from  the 
static  pressure  distribution  is  shown  on  each  of  the  heat  transfer 
graphs  for  It  is  seen  that  the  length  of  the  heat  transfer 

interaction  region  (from  the  station  at  >diich  h-r  first  rises 
above  the  zero  pressure  gradient  value  to  where  peaks)  is 
consistently  less  than  (the  pressure  interaction  length).  In 
general}  increases  linearly  to  the  peak  and  then  decays  until 
the  end  of  the  heat  transfer  instrumentation  is  reached. 

The  effect  of  various  blowing  rates  on  the  distribution  of  /r 
for  •  12°  is  shown  in  Figure  33.  Increased  blowing  appears  to 
slightly  decrease  hr  at  the  peak.  The  same  comment  applies  to 
the  ■  6®  data}  but  essentially  no  effect  is  found  at  ^  ■  9.5®. 
However}  since  the  influence  of  blowing  is  small  any  trend  will 
be  partially  masked  by  the  natural  scatter  in  the  data. 


U.  Recovery  Factor 

In  this  report  the  usual  definition  of  recovery  factor,  ^  ,  is 
enployed : 


T, 

T^-77 


where  Zw  is  the  adiabatic  wall  temperature,  %  is  the  stagnation 
te]tg)erature  of  the  free-stream,  and  T,  is  the  static  temperature 
outside  the  boundary  layer  on  the  test  plate  with  no  impinging 
shock.  For  adiabatic  flow  external  to  the  boundary  layer,  ^  , 
can  be  written  as 


( t*  O.Z 


Tr 


-  / 


O.Z  /f/ 


where  ■  2.98  for  the  nominal  Mach  3  data.  The  distribution  of 
^  along  the  test  plate  for  various  blowing  rates  and  with  ■  0 
is  shown  in  Figure  3U>  It  is  seen  that  boundary  layer  blowing  has  a 
small,-  but  explicit,  effect  on  71  with  increasing  ,  yielding 

lower  values  of  recovery  factor.  A  definite  "wiggle"  in  the 
curves  is  found  downstream  of  station  12.0  which  implies  that  a 
disturbance  exists  which  influences  the  boundary  layer. 
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In  nan/  engineering  applications  %  for  turbulent  flow  is  considered 
to  be  equal  to  the  cube  root  of  the  Prandtl  numberi  .  7(1/  is 
about  $2$9R  for  the  data  of  Figure  3ht  and  the  corresponding  Prandtl 

nunber  is  0.71j  therefore,  ^  ■  0,89U.  In  Figure  3U  this 

Talue  of  ^  is  a  reasonable  approximation  to  the  experimental 
reeover7  factor,  espeolall7  for  the  lower  blowing  rates. 

The  adiabatic  wall  temperatures  were  obtained  with  the  power  to  the 
heater  elements  turned  off  and  with  the  tunnel  operating  continuously 
for  more  than  one-half  hour.  Thus,  it  is  felt  that  steady  state 
temperatures  had  been  reached  and  the  thermocouple  outputs  were 
truly  indicative  of  the  correct  Unfortunately,  no  runs  of  this 

nature  were  made  for  •<1>  0°.  If  it  is  desired  to  gain  further 
information  regarding  ^  ,  then  it  must  be  done  indirectly  by 

investigating  the  effective  recovery  factor,  TJx  t  distribution 
which  is  defined  by  the  relationship. 


The  effective  recovery  temperature,  ,  was  defined  in  Section  II  B 
with  the  aid  of  Figure  16,  Since  the  value  of  at  a  given 
station  denotes  the  Isotheivaal  wall  temperature  tAich  corresponds 
to  sero  heat  transfer,  it  is  a  more  significant  quantity  in 
calculating  heat  flow  than  Apparently,  the  tacit  assumption 

has  been  made  by  most  investigators  that  "ti  »  although  they 
are  not  necessarily  the  same  as  pointed  out  by  Gadd  in  Reference  (8). 
However,  a  large  difference  between  the  two  temperatures  (and  hence 
between  ^  and  )  is  not  to  be  expected.  In  fact,  Gadd  reports 
that  Thomann  in  Sweden  found  and  %,  to  be  virtually  the  same 
for  experiments  performed  at  Mach  1.8*.  Distributions  of  %  on 
the  test  plate  with  the  generator  plate  removed  and  with  the 
generator  at  zero  angle  of  attack  are  shown  in  Figures  3?  and  36. 

A  comparison  between  ^  and  for  <,  ■  0  is  illustrated  in 

Figure  37.  For  •  1,0,  ^  and  are  practically  identical 


^t  should  be  noted  a  larger  experimental  error  is  to  be  expected  in 
71  than  inTi^  .71  ^  particular  station  is  the  result  of 

extrapolating  a  straight  line  determined  by  four  experimental  sets  of 
values  of  heat  transfer  and  tenperature.  Whereas,  TT./  at  the  same 
station  is  obtained  directly  from  a  single  thermocouple  reading. 
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except  in  the  region  between  stations  12.0  and  lli.O  where  a  difference 
as  large  as  0.03  is  found.  The  data  in  Figure  37  obtained  with 
bounder/  layer  blowing  indicate  a  fairly  consistent  but  small  difference 
between  \  and  has  the  sane  general  distribution  and 

is  about  0.01  less  than  \  .  Thus,  it  is  seen  that  although  \  and 
are  not  identical  they  both  illustrate  similar  distzdbutlons  along 
the  test  plate.  Therefore,  any  qualitative  trends  found  true  for 
should  also  be  valid  for  \  . 

Ihe  effective  recoveiy  factor  distributions  on  the  test  plate  for 
generator  plate  deflections  of  6.0*’,  and  12.0^  are  shown  in 
Figures  38  through  1(0.  From  all  three  figures  it  nay  be  seen 
that  the  influence  of  blowing  diminishes  with  distance  such  that 
beyond  station  12.0  the  differences  are  within  experimental  accuracy 
(about  4.008).  This  effect  is  not  due  to  the  impinging  shock  since  • 
the  data  with  the  shock  generator  plate  removed  (Figure  3^)  indicate 
the  same  general  trend.  There  is  a  small,  consistent  increase  in  ^ 
downstream  of  station  12.0  for  increasing  shock  strength  from 

0.88  for  ■  0®  to  ^s^.91  for  ^  ■  12.0°.  The  low  values 
of  upstream  of  station  12.0  for  ^  ■  9.5°  (see  Figure  39) 
apparently  resulted  from  recording  data  before  the  heat  flow  had 
reached  equilibrium  conditions  in  the  insulating  plastic  laminate. 

The  effect  of  not  quite  obtaining  steady  state  heat  flow  will 
produce  only  a  small  error  in  the  heat  transfer  and  recovery 
temperature  measurements  in  the  interaction  region,  i.e.,  downstream 
of  station  12.0. 


B.  MACH  5.0  DATA 

A  free-strean  Mach  number  of  5*02  was  determined  from  the  settling 
chamber  pressure,  ^  ,  and  from  the  pitot  tube  placed  below  the  test 
plate.  However,  the  static  pressure,  ,  measured  on  the  wind 
tunnel  wall  with  a  tap  placed  below  and  slightly  forward  of  the  test 
plate  indicated  values  about  30^  higher  than  ^  (the  static 
pressxire  calculated  by  allowing  an  isentropic  expansion  from  ^ 
until  a  Mach  number  of  5*02  is  attained).  The  static  pressure  on 
the  test  plate,  ,  was  slightly  higher  than  .  An  unrealistic 
Mach  number  of  about  li»$  was  obtained  from  the  Rayleigh  pitot  formula 
by  using  the  values  of  and  the  pitot  tube  pressure.  A  Mach 
number  slightly  greater  than  l(.8  was  calculated  by  assuming  an 
isentropic  expansion  from  the  stagnation  pressure  to  a  free- 
stream  static  pressure  equal  to  .  Of  course,  one  possible 
reason  for  the  lack  of  correlation  of  the  various  calculated  Mach 
numbers  is  a  pressure  reading  error.  Since  ^  and  vere 


17 


NORTH  AMERICAN  AVIATION  .  INC. 


COlUMlUt  OlVlilON 
COlUMIUf  U.  OHIO 


HA  62fl-795 


fairly  close  in  Talcs  it  is  reasonable  to  suspect  the  accuracy  of 
the  pitot  tube.  Rowevor,  another  explanation  may  be  based  on  the 
possible  angularity  of  the  flow  (apparently,  the  Mach  5  nozzle  was 
never  checked  for  flow  angularity).  That  is,  the  nozzle  flow 
possibly  had  not  completed  its  e:^nsion  before  reaching  the  leading 
edge  of  the  test  plate,  and  thus,  the  flow  was  con^ressed  by  the 
top  of  the  test  plate,  producing  the  increase  of  P,  over  ^  . 

If  it  is  assumed  that  the  flow  was  symmetrical,  then  the  generator 
plate  was  also  subjected  to  an  under  expanded  flow  field,  and  as  a 
result,  the  generator  surface  facing  the  test  plate  compressed  the 
flow.  There  are  several  reasons  for  accepting  this  latter  explanation. 
First,  it  was  found  that  with  the  generator  plate  at  zero  angle-of- 
attack  (i.e.,  parallel  to  the  tunnel  centerline)  there  was  a  positive 
pressure  rise  on  the  test  plate  between  model  stations  10.0  and  11.0 
(see  Figures  Ul  and  1^2).  This  is  the  region  where  a  weak  shock 
emanating  from  the  leading  edge  of  the  generator  plate  would  inplnge 
on  the  test  plate.  Second,  it  is  generally  accepted  that  the  pitot 
tube  pressure  is  insensitive  to  small  changes  in  angle-of -attack; 
consequently,  the  pitot  tube  beneath  the  test  plate  would  have 
measured  the  correct  local  Mach  number  in  a  slightly  expanding  flow 
field  if  the  pitot  pressure  was  determined  accurately.  Therefore, 
the  reason  for  the  discrepancy  between  the  measured  values  of 
and  f}  as  compared  with  the  value  of  P  determined  indirectly 
from  a  pitot  tube  measurement  is  either  that  the  pitot  tube  pressure 
was  inaccurate  or  that  the  flow  field  was  expanding  in  the  region 
of  the  test  plate  and  generator  plate  leading  edges,  thereby 
creating  an  effective  positive  angle-of -attack  (about  2®)  for  both 
surfaces.  Fortunately,  both  explanations  yield  essentially  the 
same  Mach  number  (Mi  «  U.77)  on  the  upstream  portion  of  the  test 
plate  and  very  nearly  the  same  value  of  the  calculated  reflected 
shock  invlscid  pressure  rise  on  the  test  plate  for  a  given  angle - 
of-attack.  For  the  sake  of  simplicity  all  calculations  for  the  nominal 
Mach  5  data  assume  that  the  flow  upstream  of  the  test  and  generator 
plates  is  parallel  to  the  tunnel  centerline  and  at  a  Mach  number  of 
lt.77.  However,  the  titles  of  the  graphs  refer  to  a  free-stream  Mach 
number  of  5.02  since  this  is  consistent  with  the  tabulated  data 
obtained  from  the  Pt.  Mugu  data  reduction  facilities. 


1.  Static  I^ressures  on  the  Test  Plate 

The  Mach  5  static  pressure  distributions  are  shown  in  Figures  Ul, 
U2  and  U3  for  the  various  combinations  of  <.  and  ^/p,  tested  in 

the  wind  tunnel.  The  increase  in  through  the  interaction 

is  generally  less  than  the  inviscid  value  shown  on  the  graphs  for 
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a  simple  reflected  shock.  Part  of  the  difference^  as  noted  previously, 
may  be  due  to  the  Inviscid  pressure  increase  not  being  quite  correct 
because  of  the  uncertainty  in  determining  the  true  free-stream 
conditions.  However,  for  oC  >  11°  and  12°,  this  uncertainty  cannot 
completely  explain  why  a  greater  inviscid  pressure  rise  was  attained. 
The  actual  reason  is  that  for  the  two  highest  values  of  od  the 
expansicm  emanating  from  'Uie  generator  plate  trailing  edge  impinged 
on  the  test  plate  in  a  region  of  increasing  pressure  and  did  not 
allow  any  further  increase.  Instead,  the  pressure  decayed  quite 
rapidly.  From  the  results  of  Bogdonoff,  Reference  (2),  it  is  believed 
that  the  boundary  layer  on  the  test  plate  upstream  of  the  expansion 
iiiq>ingement  point  was  not  influenced  by  the  expansion. 

During  all  Mach  5  tests  a  wide  strip  of  .033"  grit  was  placed  1.5" 
fnom  the  leading  edge  of  the  test  plate.  Apparently  when  no  boundary 
layer  blowing  was  applied  (note:  l.U  implies  no  blowing 

since  1»U),  the  flow  was  laminar  despite  the  effect  of  the 

grit.  The  pressure  distributions  with  no  blowing  for  •  10° 
and  11°  indicate  a  comparatively  slow  rise  and  then  a  flattening 
out  before  the  pressure  increases  rapidly  to  its  peak  value.  This 
type  of  curve  is  indicative  of  flow  which  separated  while  laminar 
and  then  became  turbulent,  e.g.,  see  Reference  (9).  That  the  flow 
was  laminar  for  ■sc.  l.h  is  further  verified  by  the  heat  transfer 

data  discussed  later.  Air  injection  into  the  boundary  layer  produced 
transition  on  the  upstream  portion  of  the  plate  as  evidenced  by  the 
shorter  interaction  length  e>C  •  11°  and  the  increase  in  heat 
transfer  upstream  of  the  interaction.  Increased  blowing  beyond 
that  necessary. to  obtain  turbulent  flow  had  a  negligible  effect 
(within  the  accuracy  of  the  data)  on  the  interaction  length  for 
a  given  cC  . 

With  boundary  layer  blowing,  the  three  generator  angles,  *>C  ■  9.7°, 

11°,  12®,  yield  attached,  indpiently  separated,  and  separated  flow, 
respectively.  Incipient  separation  is  considered  to  occur  11® 

because  a  slight  inflection  in  the  pressure  distribution  is  found 
within  the  interaction  region.  It  is  probable  that  if  the  pressure 
taps  were  more  closely  spaced,  incipient  separation  would  have  been 
determined  at  a  few  tenths  of  a  degree  less  than  11°.  At  ■  12® 
a  distinct  Inflection  in  the  pressure  distribution  in  the  interaction 
region  Infers  that  the  flow  was  separated. 
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2.  Boundary  Layer  Rakes 

The  Mach  5  velocity  profiles  at  stations  9.0  and  15.0  determined 
from  the  pitot  tube  boundaiy  layer  rake  pressures  are  shown  In 
Figures  hh,  U5  and  U6.  With  no  boundary  layer  blowing,  i-ii. 

the  profiles  at  station  9.0  indicate  a  region  within  the  boundaxy 
layer  for  which  the  slope  of  the  profile,  ,  is  practically  zero. 

In  addition,  the  data  points  within  a  run  (l.e.,  keeping  all  conditions 
the  same  but  changing  the  wall  temperature)  gave  different  profiles, 
which  Implies  either  an  unusual  sensitivity  to  temperature  or,  more 
likely,  unsteady  flow.  With  blowing  the  shape  of  the  velocity 
profiles  appear  normal  for  «  0®  and  with  the  generator  plate 
removed.  The  profiles  at  stations  9.0  and  15.0,  however,  are  almost 
Identical,  Indicating  a  constant  thickness  of  O.3I1''.  It  would  seem 
that  the  only  straightforward  reason  for  the  boundary  layer  not 
thickening  with  distance  Is  the  result  of  a  favorable  pressure 
gradient  or  a  three-dimensional  relieving  effect.  Neither  explanation 
appears  to  be  applicable  In  this  case  since  the  pressure  distribution 
Is  flat  and  the  oil  flow  studies  of  the  streamlines  showed  a  two- 
dimensional  pattern. 

For  >  0®  the  upstream  probe  was  unaffected  by  the  interaction 
except  at  0^  ■  12°.  At  «  12®  the  interaction  started  near  station 
8.5  and  the  forward  velocity  profile  was  somewhat  distorted,  which 
Indicates  an  Increase  In  the  velocity  defect  as  compared  to  the 
profiles  for  12®.  For  OP  all  velocity  profiles  obtained 
with  the  downstream  probe  were  similar  in  shape.  The  velocity 
defect  was  less  than  the  defect  at  station  9.0  (the  converse  effect 
was  fouiK]  at  Mach  3)  and  0  between  ^  ■  0.05"  and  ^  «  0.085". 

The  value  of  the  velocity  ratio,  ^Aje.  ,  at  which  'X  0, 

Increased  with  Increasing  shock  strength.  * 

Increasing  blowing  from  li  to  9.5  had  no  consistent 

effect  on  the  shape  of  the  velocity  profiles,  and  the  boundary 
layer  thickness  remained  at  about  0.3U". 


3.  Heat  Transfer  Coefficient 

The  heat  transfer  coefficient  distributions  for  Mach  5  are  shown 
In  Figures  U7  through  60.  Two  theoretical  lines  corresponding  to 
zero  pressure  gradient  heat  transfer  coefficient  are  drawn  on  each 
graph.  One  line  was  derived  from  the  straightforward  use  of  the 
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turbulent  heat  transfer  theory  of  Reference  (7)  assuming  the  boundary 
layer  becomes  turbulent  at  the  leading  edge  of  the  test  plate.  The 
second  line  was  also  calculated  from  the  theory  of  Reference  (7),  but 
was  adjusted  to  account  for  the  experimentally  determined  momentum 
thickness.  Unlike  the  Kach  3  case,  the  momentum  thickness  was  not 
altered  significantly  by  varying  f^e/p  above  the  minimum  value 

necessary  to  obtain  turbulent  flow.  Consequently,  only  one  line  was 
needed  to  show,  the  theoretical  heat  transfer  distribution  based  oh 
measured  6  .  It  may  be  seen  that  the  strict  use  of  the  theory  yields 
better  agreement  with  the  data  than  the  theory  ad justed^ to  account  for 
measured  0  ,  ^ 

The  length  of  the  pressure  interaction  region,  ,  is  noted  on  the 
heat  transfer  graphs  for  >  0®.  Consistent  with  the  trend  found 
at  Mach  3j  the  heat  transfer  interaction  length  is  less  than  . 
Also,  in  agreement  with  the  Mach  3  results,  an  increase  in  blowing 
rate  will  either  have  no  effect  or  will  increase  the  heat  transfer 
interaction  length  and  slightly  decrease  the  peak  heat  transfer 
coefficient.  Any  definite  trend  is  hidden  by  scatter  In  the  data. 
Upstream  of  the  Interaction  with  no  boundary  layer  blowing  (see 
Figures  $8  and  $9)  the  value  of  Ar  is  quite  low  which  Is 
indicative  of  laminar  or  the  initial  phase  of  transitional  flow. 
Apparently  transition  to  turbulence  occurs  rapidly  In  the  interaction 
region  as  implied  by  the  hj-  distribution. 


l(.  Recovery  Factor 

The  Mach  $  adiabatic  wall  temperature  distribution  on  the  test  plate 
was  determined  for  cC  “  7.U®,  10°  and  12°  by  recording  the 

thermocouple  outputs  after  the  tunnel  had  been  operating  for  at 
least  30  minutes.  The  results  In  terras  of  recovery  factor,  71  , 
are  shown  in  Figures  61,  62  and  63.  In  Figure  61  It  is  seen  that 
for  -  0®  and  no  boundary  layer  blowing  ^  increases  slowly  with 
increasing  distance  from  the  leading  edge  and  then  peaks  and  begins 
to  decrease  near  station  10.5.  A  similar  variation  of  7Z  with 
distance  is  found  in  the  data  obtained  by  Brlnich  (see  Reference  10) 
with  a  hollow  blunted  cylinder  at  Mach  5.  Brlnich  showed  that  for 
a  transitional  boundary  layer  ^  goes  through  a  peak  with  the  peak 
in  corresponding  to  the  end  of  transition  (this  was  verified 
with  Schlieren  photographs).  Bluntness  also  raised  the  entire 
level  of  the  ^  distribution  and  delayed  transition  when  compared 
with  the  results  obtained  with  a  sharp  leading  edge  cylinder.  The 
value  of  the  peak  in  Figure  6l  is  0.92  which  compares  favorably  with 
^  ■  .912  found  by  Brlnich  when  the  semi-width  of  the  bluntness 
was  .030".  (It  is  recalled  that  the  height  of  the  grit  employed  in 
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the  present  test  was  .033".)  Also,  the  rate  of  change  of  ^  with 
distance  agrees  well  with  the  hollow  cylinder  tests  up  to  the  peak, 
but  Brinich's  data  Indicate  a  slower  decay  in  V,  downstream  of  the 
Biaximun  1l  .  Thus,  in  the  present  tests  it  appears  that  without 
boundary  layer  blowing,  grit  produces  the  same  qualitative  effect  on 
the  boundary  layer  as  a  blunted  leading  edge. 

As  was  mentioned  previously,  the  combined  effects  of  grit  and  boundary 
layer  blowing  produced  tranirltion  upstream  of  the  instrumented  portion 
of  the  test  plate.  This  was  evidenced  by  the  pressure  and  heat  transfer 
data.  For  cC  •  0®,  ^  decreased  slowly  with  distance  as  would  be 
expected  by  analogy  with  the  hollow  cylinder  data  of  Reference  (10). 
Increasing  the  blowing  rate  from  ^n/f,  ■  3.9  to  •  8.9  had 

practically  no  effect  on  7^  .  *  * 

In  Figures  62  and  63  an  unexpected  distribution  of  ^  is  found  for 
eC  y  0®.  ^  suddenly  increases  and  then  decreases  within  a  space 
of  two  Inches  or  less.  Near  the  end  of  the  instrumented  region 
(between  stations  lli.O  and  16.0)  approximately  the  same  value, 
fi  ■  0.883,  Is  attained  whether  or  not  blowing  is  applied.  In  the 
same  region  0.892  for  •<-  0®.  The  qualitative  effect  of 

shockwave  boundary  layer  interaction  on  ^  at  Mach  5  Is  sketched 
in  Figure  70.  At  Mach  3  this  unusual  variation  in  ^  did  not  occurj 
however,  the  effect  was  found  at  Mach  3.1  by  Brlnlch  (Reference  11) 
when  two-dimensional  roughness  elements  were  placed  on  a  hollow 
cylinder.  As  noted  at  the  bottom  of  Figure  70,  a  wedge  type  roughness 
element  with  the  flat  surface  facing  the  oncoming  flow  produces  the 
same  type  of  variation  \xi  ?i  as  found  with  the  present  Mach  5  data 
for  oL>  0®.  Brlnlch  also  found  that  comparatively  large  variations 
in  %  occur  in  the  vicinity  of  a  roughness  element  placed  downstream 
of  one  or  more  (the  maximum  tested  was  four)  similar  elements.  That 
is,  the  effect  may  be  induced  in  the  same  boundary  layer  more  than  once. 
It  is  also  interesting  to  note  that  in  an  experiment  conducted  by 
Oadd  at  Mach  2,hh  a  result  similar  to  BrinichJs  and  the  present  Mach  5 
data  was  found.  In  Gadd's  experiment  (Reference  12)  a  slab  was  placed 
in  the  turbulent  boundary  layer  on  a  wind  tunnel  wall.  The  height 
of  the  slab  was  about  2.5  times  the  undisturbed  boundary  layer 
thickness,  ,  whereas  in  Brinich's  test  the  roughness  elements  were 
probably  less  than  one-half  S  .  Oadd  found  that  the  recovery 
factor  increased  and  peaked  upstream  of  the  slab  and  decreased  down¬ 
stream  of  the  rearward  facing  step  to  a  value  below  the  undisturbed 
boundary  layer  recovery  factor.  Why  no  large  change  in  ^  was 
observed  for  the  present  Mach  3  tests  is  not  understood,  but  it  may 
be  related  to  the  absence  of  grit.  It  would  have  been  interesting 
to  see  if  this  phenomenon  would  have  occurred  at  Mach  5  if  the  test 
plate  were  long  enough  to  allow  transition  without  the  aid  of  grit. 
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It  nay  be  that  grit  produces  a  change  in  the  turbulent  mechanism 
such  that  further  disturbances  to  the  boundary  layer  yield  a  somewhat 
different  result  than  would  be  found  for  a  boundary  layer  which 
tuiderwent  natural  transition.  This  point  is  discussed  further  with 
regard  to  the  length  of  the  interaction  region  in  Section  IV  A. 

The  effective  recovery  factor,  ,  was  defined  in  Section  III  A  li 
as  a  function  of  the  effective  recovery  temperature,  TT  •  The 
distribution  of  9^^  on  the  test  plate  for  various  values  of 

and  •(  is  shown  in  Figures  6U  through  69.  Unfortunately,  the  method 
of  obtaining  TJ  (see  Figure  16)  was  rather  Inaccurate,  especially 
in  regions  of  low  heat  transfer.  Thus,  for  ■  0®  or  upstream  of  ^ 
the  interaction  for  >  0®  a  large  amount  of  scatter  is  expected 
in  .  Inspection  of  Figure  6ii  reveals  that  considerable  scatter 

does  exist  in  ^  for  the  undisturbed  boundary  layer  (about  +.02). 

In  the  region  of  comparatively  high  heat  transfer  more  accurate 
results  are  obtained.  In  Figures  65  through  68  this  is  evidenced 
by  the  data  downstream  of  station  11.0.  The  relationship  between 
^  and  as  Illustrated  in  Figure  69,  is  the  same  as  that  found 

at  Mach  3»  i.e.,  yields  a  distribution  similar  to  but  the 
variation  is  between  wider  limits. 


( 
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17.  DISCUSSION  AND  CORREUTION  OF  DATA 


A.  EFFECT  OF  BOUNDARY  LAIER  TRIPS  ON  CORRELATION  OF  PRESSURE 
INTERACTION  LENGTH 

In  addition  to  the  present  data,  pressure  data  from  References (1), 

(2)  and  (6)  were  used  In  the  correlation  of  the  pressure  Interaction 
length,  .  (The  applicable  pressure  data  from  Reference  (1)  are 
presented  in  Figures  73  and  7U.)  In  the  present  tests,  the  pressure 
taps  and  the  boundaiy  layer  rakes  were  not  located  on  the  heated 
portion  of  the  test  plate.  The  pressure  data  from  the  foregoing 
references  were  from  tests  employing  adiabatic  walls.  As  shown  in 
References  (8)  and  (13)j  however,  flow  properties  are  fairly  insensitive 
to  small  differences  between  7^  and  »  Therefore  it  is  assumed 

that  small  rates  of  heat  transfer  have  negligible  effects  on  , 

&!  *  ^*/p,  and  A'W  .  Also,  since  the  heat  transfer  elements  in 
this  test  were  located  in  a  small  section  centered  on  the  test 
plate,  it  is  assumed  that  the  heat  transfer  did  not  introduce 
significant  three-dimensional  effects  into  the  flow. 

Two  correlations  ot  A 'ii  were  made  and  are  shown  in  Figures  71  and 
72.  The  data  used  in  making  these  plots  are  listed  in  Table  II.  As 
noted  previously,  the  boundary  layer  thicknesses  were  obtained  from 
velocity  profiles  by  extrapolation  of  the  straight  line  portion  of 
the  log- log  plots  of  ^  versus  through  •  1.  These  values 

of  S  were  then  adjusted  to  account  for  the  length  of  boundary  layer 

run  from  the  point  where  the  measurements  were  taken  to  the  point 
where  the  pressure  interaction  began.  To  make  this  adjustment  it 

was  assumed  that  S  varied  as  %  '  where  is  measured  from 
the  leading  edge  of  the  plate.  Since  the  flow  did  not  become 
turbulent  at  the  leading  edge  in  this  or  the  experiments  of 
References  (1)  and  (6),  this  adjustment  was  somewhat  in  error  even 
if  the  ^  dependence  assximed  was  accurate.  The  momentum  thicknesses 
were  obtained  from  the  probe  measurements  and  adjusted  for  distance 
as  previously  described  using  the  same  ■oc  dependence  to  the  beginning 
of  the  interaction.  The  pressure  interaction  length,  A'K  ,  was 
defined  to  be  the  distance  from  the  station  where  the  pressure  first 
began  to  rise  above  P,  to  the  pressure  peak,  .  The  quantity 

may  be  used  to  characterize  the  pressure  interaction. 
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The  variations  of  and  with  this  parameter  are  shown 

In  Figures  71  and  72.  It  is  seen  that  and  vary  almost 

linearly  with  the  interaction  parameter. 

Figure  71  shows  two  pressure  interaction  correlation  lines.  On  the 
lower  line  lie  the  present  •  2.98  results  and  the  experimental 
data  of  Bogdonoff  (Reference  2)  and  Simon  (Reference  1).  On  the 
upper  line  are  found  the  present  ■  U.77  data  and  Kuehn's  data 
(Reference  6).  From  Table  II  it  may  be  seen  that  in  comparing 
the  present  Mi  ■  li.77  data  and  Kuehn’s  data,  and  Mi  differ  by 
about  a  f. actor  of  2  and  1.6,  respectively,  and  although  these  two 
sets  of  data  correlate  well,  they  do  not  correlate  with  the  remaining 
data  which  define  the  lower  line  in  Figure  71.  Since  the  data  on 
the  lower  line  were  obtained  from  experiments  covering  a  range  of 
values  of  ^  and  Mi  it  appears  that  the  separation  of  all  data 
points  into  two  separate  correlation  lines  must  be  attributed  to 
the  method  of  inducing  transition.  For  the  Mi  ■  U.77  data  of  this 
test  and  Kuehn's  experiments,  leading  edge  boundary  layer  trips 
were  used.  In  the  present  case  transition  was  affected  by  a 
combination  of  grit  plus  boundary  layer  blowing,  whereas  Kuehn  forced 
transition  by  placing  a  rearward  facing  step  near  the  leading  edge  of 
an  otherwise  flat  plate.  It  should  be  noted  that  Bogdonoff 's 
measurements  were  conducted  along  a  wind  tunnel  wall  where  natural 
transition  occurred.  Figure  71  shows  that  the  two  methods  of 
artificially  tripping  the  boundary  layer  have  the  same  effect  on 
»  however,  a  gross  difference  is  indicated  between  the 

cases  where  transition  is  induced  as  compared  with  the  flows  that 
experience  natural  transition. 

The  second  interaction  length  correlation.  Figure  72,  represents 
a  plot  of  versus  It  is  seen  that  in  this 

case  three  distinct  lines  are  determined  by  the  data.  Employing 
the  reasoning  given  above,  it  is  plausible  to  attribute  the  grouping 
of  the  points  to  the  method  of  cr.  hieing  transition.  It  is  not  too 
surprising  that  the  parameter  is  more  sensitive  to  the 

details  of  the  boundary  layer  :  c  than  ,  since  Q 

is  certainly  more  sensitive  thii".  to  the  past  history  of  the 
boundary  layer  and  any  significant  changes  in  the  turbulent 
mechanism.  It  should  be  noted  that  attempts  to  correlate  the 
interaction  length  using  parameters  involving  local  skin  friction 
upstream  of  the  Interaction  yielded  the  same  type  of  results  as  those 
discussed  previously.  A  correlation  which  is  useful  in  predicting 
interaction  length  for  a  given  set  of  upstream  conditions  is 
considered  in  Section  V. 
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The  most  significant  fact  to  be  gained  from  the  foregoing  discussion 
is  that  the  method  of  inducing  transition  may  have  a  noticeable  effect 
on  the  subsequent  behavior  of  the  turbulent  boundary  layer.  This  was 
evidenced  by  the  two  correlations  of  shown  in  Figures  71  and  72. 
In  addition,  it  is  also  important  to  note  that  increasing  the  boundary 
leyer  blowing  ratio,  ^/pt  ,  above  the  minimum  value  needed  to  induce 
turbulence  does  not  alter  the  correlation.  That  is,  a  particular 
correlation  cannot  differentiate  between  changes  in  S,  or 
produced  by  (1)  varying  the  free-stream  conditions  (including  the 
length  of  boundary  layer  runjor  (2)  merely  varying  the  rate  of  air 
injection  in  the  boundary  layer.  Thus,  it  appears  (although  it 
certainly  has  not  been  proven)  that  if  the  boundary  layer  can  be 
thickened  with  blowing  without  producing  a  distorted  velocity  profile, 
then  it  will  interact  with  a  shock  wave  in  exactly  the  same  manner  as 
a  boundary  layer  thickened  the  same  amount  by  either  decreasing  the 
Reynolds  number  per  unit  length  or  increasing  the  length  of  run. 

This  irqjlies  that  leading  edge  blowing  may  be  used  to  simulate 
Reynolds  number  effects. 


B.  STANTON  NUMBER  RATIO  DISTRIBUTION 

In  order  to  correlate  the  heat  transfer  data  of  the  present  test  in 
terms  of  non-dimensional  quantities,  one  reasonable  parameter  to 
choose  is  the  local  Stanton  number. 


where  is  the  local  mass  flow  external  to  the  bo\mdary  layer. 

Increasing  shock  strength  (i.o.,  oA  going  from  0°  to  12®)  infers 
higher  values  of  in  the  shockvrave  boundary  layer  interaction 

region.  However,  if  the  viscous  flow  maintains  itself  as  a  zero 
pressure  gradient  boundary  layer  and  merely  adjusts  to  the  change 
in  the  external  inviscid  flow  produced  by  the  reflected  shock,  then 
despite  the  increase  in  would  also  increase  ^ut  only 

slightly).  The  increase  is  primarily  due  to  the  change  in  Mach 
number  through  a  shock  and  the  attendant  effect  on  the  zero  pressure 
gradient  boundary  layer  (see  Reference  22).  Therefore,  an  appreciable 
change  in  in  the  interaction  region  is  an  indication  of  a 

corresponding  change  in  tJie  boundary  layer.  To  determine 
certain  assumptions  must  be  made  in  the  calculations . 
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The  flow  properties  external  to  the  boundary  layer  are  known  only 
indirectly  from  the  static  pressure  distribution  on  the  test  plate. 

The  invariance  of  static  pressure  in  the  ^  direction  through  a 
boundary  layer,  derived  from  an  order  of  magnitude  argument,  is 
no  longer  valid  in  the  interaction  region.  The  large  streaniwise 
gradients  of  static  pressure  on  the  wall  are  associated  with 
streamwise  curvature  in  the  flow  which,  in  turn,  infers  the 
existence  of  vertical  pressure  gradients.  Thus,  there  will  be  some 
discrepancy  between  the  test  plate  static  pressure,  P  ,  and  the 
static  pressure,  ^  ,  in  the  inviscid  flow  external  to  the  boundary 

layer.  Nevertheless,  P  is  a  reasonable  approximation  to  /2  and  no 
significant  differences  of  a  qualitative  nature  are  to  be  expected. 

An  additional  assumption  is  made  that  at  each  point  in  the  inviscid 
external  flow,  the  local  flow  properties  Including  the  static  pressure 
are  the  result  of  the  undisturbed  flow  on  the  test  plate  going 
through  a  reflected,  oblique,  inviscid  shock-M'.  In  Figures  7$  and  76 
Mach  number,  mass  flow,  ratio,  and  initial  flow  deflection  angle 
are  plotted  versus  pressure  ratio  across  a  reflected  oblique  shock- 
wave  in  inviscid  flow  for  upstream  Mach  numbers  of  2.98  and  U.77, 

These  two  graphs  were  employed  in  the  calculation  of  Stt  from  the 

heat  transfer  and  pressure  data.  For  each  experimental  value  of 
pressure  ratio  a  mass  flow  ratio  was  determined  from  Figure  75  or  76. 
The  local  Stanton  number  ratio  was  calculated  using  this  mass  flow 
ratio  and  the  local  experimental  value  of  4r  .  The  results  of  the 
calculations  are  given  in  Figures  77  through  80.  In  these  figures 
is  plotted  versusf'’^-‘>'i^-5^  where  the  zero  pressure 

gradient  Stanton  number  St,  ,  the  pressure  interaction  length “5^  , 

and  the  station  at  which  the  interaction  begins  are  also 

experimentally  determined.  Each  point  in  the  graphs  corresponds 
to  the  location  of  a  static  pressure  tap  or  the  center  of  a  heat 
transfer  element.  Since  in  no  case  are  the  two  locations  coincident 
it  was  nece.ssaiy  to  fair  the  pressure  data  to  obtain  a  pressure  at  the 
center  of  a  hest  transfer  element  and  conversely  to  obtain  the  heat 
transfer  coefficient  at  a  pressure  tap  location.  In  all  cases  the 
fairings  were  made  through  the  actual  data  points  even  when  the  data 
appeared  somewhat  erratic.  However,  a  certain  amount  of  arbitrariness 
remained  while  fairing  in  a  region  where  large  pressure  or  heat 
transfer  gradients  occurred.  It  is  seen- that  a  certain  rough 
similarity  exists  among  all  four  graphs.  The  Stanton  number  ratio  is 


A 

This  assumption,  of  course,  is  never  quite  true  except  at  the  end  of 
the  pressure  interaction.  However,  since  the  flow  external  to  the 
boundary  layer  will  undoubtedly  experience  an  increase  in  entropy,  even 
near  the  beginning  of  the  interaction,  it  is  probably  more  valid  to 
relate  the  measured  static  pressures  to  a  reflected  shock  system  than 
to  use  the  sinqjle  isentroplc  relationships  (as  Qadd  did  in  Reference  8). 
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less  .than  unity  in  the  upstream  portion  of' the  interaction  region. 

It  rises  steeply  and  peaks  1  ,  and  then  decays  at 

a  comparatively  slow  rate.  The  lowest  value  of  is  less 

than  O.U  and  the  peak  is  between  1.5  and  1.8.  The  overall  variation 
In  •  1-8  given  in  Figure  8l.  Certainly,  the  area  covered  by 

the  shading  in  Figure  8l  is  too  great  to  claim  that  the  particular 


set  of  parameters 


and 


chosen  uniquely  describe 


the  heat  transfer  phenomenon  downstream  of  ,  However,  there  is 
a  considerable  amount  of  unavoidable  error  in  determining  the  correlation 
parameters  which  may  be  responsible  for  almost  all  of  the  spread  in 
the  data.  At  Mo  »  2.95  the  Inability  to  accurately  determine 'Xi 
and  40S  (due  to  the  spacing  of  the  pressure  taps)  and  at  M.  ■  5*02 
the  errors  associated  in  calculating  (primarily  due  to  the  low 
heat  transfer  rates  caused  by  low  tunnel  densities)  are  the  most 
probable  sources  of  scatter. 


Since  a  literature  survey  did  not  reveal  any  turbulent  boundary  layer 
heat  transfer  investigations  involving  an  impinging  oblique  shock, 
no  direct  comparison  with  other  experiments  can  be  made.  However,  it 
is  of  interest  to  see  if  the  heat  transfer  associated  with  a  turbulent 
boundary  layer  disturbed  by  means  other  than  an  impinging  shock  can  be 
related  to  the  present  data.  In  Figure  82  two-dimensional  turbulent 
flow  heat  transfer  distributions  obtained  from  three  separate  tests 
are  shown  in  terms  of  Stanton  number  ratio.  A  few  comments  concerning 
the  data  in  the  figure  are  warranted.  In  certain  respects,  Oadd's 
data  (References  8  and  12),  which  were  obtained  by  disturbing  the 
flow  on  a  wind  tunnel  wall  with  a  rectangular  slab,  display  a 
similarity  to  the  present  results.  Upstream  of  the  slab, 
decreases  to  a  value  less  than  imity  in  the  separated  region  and  then 
Increases  sharply.  Again,  downstream  of  the  slab  it  is  found  that 
<  1.0  and  as  the  flow  reattaches  to  the  wind  tunnel  wall, 

rapidly  approaches  a  peak  value  of  1.7  before  It  begins  to 

decay.  Seban's  low  speed  data  (Reference  18)  in  Figure  82  show  the 
result  of  disturbing  a  flat  plate  flow  with  an  1/8"  diameter  rod 
placed  spanwlse  on  a  plate.  It  is  observed  that  the  upstream 
Influence  of  the  rod  produced  a  small  oscillation  in  the 
curve  which  is  also  fouixi  in  Gadd's  data.  The  peak  value  of 

is  approximately  1.5  and  downstream  decay  from  this  value  is  quite 
slow,  decreasing  to  a  value  of  1.2  at  a  distance  of  95  rod  diameters 
downstream  of  the  rod.  (Note  that  in  the  case  of  the  rod,  the  quantity 
//  ,  which  non-dlmenslonalizes  distance,  was  chosen  arbitrarily  to  be 
0.60"  so  that  the  data  would  fit  on  the  graph.)  The  heat  transfer  data 
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obtained  by  Naysmlth  (Reference  15)  on  a  wind  tunnel  wall  downstream 
of  a  ramp  indicates  qualitative  but  not  quantitative  agreement  with 
the  data  discussed  previously.  The  distribution  from 

Reference  (15)  in  Figure  82  does  not  contain  the  data  taken  upstream 
of  the  boundary  layer  reattachment  point  because  the  measurements  in 
the  separated  region  directly  behind  the  ramp  were  reported  to  be 
unreliable.  It  appears  that  a  peak  value  of  >3.0  may  occur 

in  the  flow  reattachment  zone.  This  maximum  value  of  is 

considerably  higher  than  what  was  found  in  the  other  Stanton  number 
ratio  distributions  (Figures  8l  and  82).  It  is  difficult  to  determine 
if  the  high  values  of  found  by  Naysmith  are  to  be  expected 

since  no  other  heat  transfer  experiments  involving  shockwave  turbulent 
boundary  layer  interaction  on  a  flat  plate  have  been  discovered  in 
the  literature.  However,  use  may  be  made  of  body  of  revolution  data 
under  certain  conditions.  If  the  Mach  number  is  fairly  high  and  the 
boundary  layer  dimensions  are  small  compared  with  the  local  body 
diameter,  the  flow  field  in  the  vicinity  of  a  body  of  revolution  is 
approximately  two-dimensional  (in  the  planar  sense).  This  is 
especially  true  at  the  juncture  of  a  cylindrical  mid-body  and  a 
conical  flare.  Two  suitable  heat  transfer  experiments  conducted  in 
a  wind  tunnel  on  pointed -nosed  bodies  of  revolution  are  reported  in 
References  (16)  and  (17).  In  Reference  (16)  the  tests  were  performed 
at  M,  ■6.8  on  bodies  with  10®  and  30®  conical  flares.  Natural 
transition  of  the  boundary  layer  occurred  on  the  cylinder  upstream 
of  the  flare  for  the  high  Reynolds  number  runs.  Peak  values  of 

reported  in  Reference  (16)  are  3.5  and  1.25*  for  the  30° 

and  10®  conical  flares,  respectively.  In  Reference  (17)  turbulent 
heat  transfer  tests  were  performed  at  Mo  ■  li.98  on  bodies  with  10° 
and  2U®  conical  flares.  A  boundary  layer  trip  on  the  nose  induced 
turbulent  flow  well  upstream  of  the  influence  of  the  cylinder  flare 
juncture.  The  peak  values  of  are  approximately  2.2  and 

l.U  for  the  21;®  and  10®  flares.  It  is  to  be  noted  that  in  both 
References  (l6)  and  (17)  the  turbulent  data  indicated  little  or  no 
separation  in  the  vicinity  of  the  cylinder  flare  juncture. 


*For  the  10®  flare  the  peak  value  of  was  found  to  be 

approximately  1.35  if  the  measured  pressures  are  employed  in 
determining  the  local  mass  flow,  ,  external  to  the  boundary 

layer.  In  Reference  (16)  was  determined  assuming  the  flow 

external  to  the  boundary  layer  identical  to  inviscid  planar  flow 
deflected  through  a  wedge  angle  equal  to  the  flare  angle. 
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Peak  values  of  ^/st,  >  have  been  found  in  three  different 
eases.  Whereas,  the  remaining  data  gave  =  1.8  (note:  the 

subscript  "2"  applied  to  Stanton  number  denotes  the  peak  value). 

It  was  noticed  that  1.8  occurred  whenever  -the  mass  flow 

ratio  in  the  flow  external  to  the  boundary  layer  was  approximately 
invariant  with  pressure  ratio,  i.e.,  whenever 

cL  (  %) 

A  graph  of  the  mass  flow  ratio  versus  the  pressure  ratio  across  a 

reflected  shock  at  various  Mach  numbers  is  shown  in  Figure  83. 

Similar  shaped  curves  are  also  determined  if  a  single  oblique  shock 

is  considered  as,  for  example,  at  the  juncture  of  a  cylinder  and  a 

conical  flare.  In  the  figure  is  defined  as  the  value  of  pressure 

v  n/B 

ratio  at  the  "elbow"  of  the  mass  flow  ratio  versus  pressure  ratio 
curve.  Where  the  elbow  actually  occurs  is  somewhat  arbitrary. 

In  all  cases  mentioned  previously  for  which  —  1.8  in  a 

supersonic  flow  field,  the  pressure  ratio  at  the  peak  Stanton  number 
location  is  less  than  foregoing  considerations 

a  tentative  expression  for  the  peak  Stanton  number  ratio  may  be 
given  as 


St.  _  ^  P  (. 

‘  '  5  (H) 

where  ^  P  (  ^  ^  ^*^)i  end  the  tacit 

assumption  is  made  that  the  inviscid  flow  properties  at  the  end  of 
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the  pressure  Interaction  (subscript  "3")  are  approximately  equal 
to  those  at  the  point  where  the  Stanton  number  peaks  (subscript  "2”). 
Flrom  the  sketch  below  It  is  seen  that  P(^j/p)  expressed  as 


where 


l.Ii<  Ki<C  (except  for  values  of  ~  1»0)  and 

;  0.lU5«  The  expression  for  was  obtained  from  the 


aerlmental  Peak  Stanton  Number  Ratio  Versus  a  Theoretical 


Inviscid  Pressure  Ratio  Increment 


Pressure  Ratio  Difference  at  Test  Mach  Number  •~' 
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a88uiq>tion  that  for  (note:  Is 

appiroxlinatel7  constant  tinder  this  condition)  the  peak  Stanton  nnnber 
ratio  varies  linearly  with  pressure  ratio.  The  constant  Ka.  vas 
determined  primarily  from  the  30°  conical  flare  data  of  Reference  (16). 
The  foregoing  equations  tor^/rt,  are  purely  speculative,  hut  it  is 

felt  that  they  represent  the  most  sinqple,  and  yet  physically  reasonable, 
formulation  that  could  be  derived  from  the  available  data.  In  order 
to  obtain  more  accurate  expressions  additional  heat  transfer  tests 
are  required. 

The  problem  of  relating  the  increase  in  Stanton  number  ratio  within 
the  shockwave  boundary  layer  interaction  region  to  a  corresponding 
change  in  the  boundary  layer  heat  transfer  mechanism  utilizing  the 
available  information  can  only  be  attacked  indirectly.  It  has  been 
conjectured  by  Naysmlth  (Reference  1$)  and  others  that  a  reattaching 
separated  boundary  layer  acts  as  if  it  contained  a  thin  "inner  boundary 
layer"  idiich  is  initiated  at  the  reattachment  point.  This  concept  is 
in  qualitative  agreement  with  ejq^riments  (References  3,  15  and  18) 
in  which  high  heat  transfer  rates  (supposedly  due  to  the  thinness 
of  the  "inner  boundary  layer")  and  low  shearing  stress  at  the 
wall  (supposedly  due  to  the  low  velocity  external  to  the  "inner 
boundary  layer")  were  found.  Further  verification  of  the  "inner 
boundary  layer"  notion  for  unseparated  flow  is  found  in  References  (16) 
and  (17)  where  experimental  supersonic  turbulent  heat  transfer  data 
on  axisymmetric  bodies  with  conical  flare  afterbodies  are  presented. 

In  both  references  it  was  found  that  reasonably  good  agreement  of 
theory  with  data  was  obtained  on  the  flare  if  the  boundary  layer  upstream 
of  the  flare  is  neglected  and  a  new  zero  pressure  gradient  boundary 
layer  beginning  at  the  cylinder-flare  Juncture  is  assumed. 

If  the  "inner  boundary  layer"  concept  is  a  valid  representation  of 
the  structure  of  the  disturbed  boundary  layer,  then  an  increase  in  fft. 

(the  recovery  factor  based  on  local  flow  conditions  external  to  the 
entire  boundary  layer)  above  the  value  found  for  zero  pressure  gradient 
flow  is  to  be  expected  downstream  of  the  beginning  of  the  "inner 
boundary  layer".  The  reason  that  an  "inner  boundary  layer"  Implies 
an  Increased  is  as  follows:  It  can  be  shown  (e.g.,  see 
References  10  and  20}  that  for  a  turbulent  boundary  layer  the  true 
local  recovery  factor,  ,  is  essentially  invariant  for  moderate 
changes  in  Reynolds  number  and  Mach  number  as  long  as  it  is  based 
on  conditions  which  are  local  frcm  the  "viewpoint"  of  the  boundary 
layer.  For  the  usual  boundaiT-  layer  •  Since  the  "inner 

boundary  layer"  is  Immersed  in  a  thicker  viscous  shear  layer  the 
true  local  Mach  number,  /*4.  ,  seen  by  the  "inner  boundary  layer"  will 
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be  less  than  the  local  Kaeh  ininber,  /f*  ,  external  to  the  entire 
boimdarjr  la/ar.  The  adiabatic  wall  temperature  for  the  case  of 
an  "inner  boondaty  layer"  is  found  from  the  relation 


tf  '  HO.iht 


If  the  fundamental  change  in  the  turbulent  viscous  flow  due  to 
its  Interaction  with  a  shockwave  is  the  creation  of  an  "inner  boundary 
layer"  (i.e.,  no-  in^rtant  changes  in  the  transport  properties),  then 
ylj,  at  a  given  point  in  the  flow  field  should  remain  unaffected 
by  the  Interaction.  Therefore,  should  increase  downstream  of 

the  interaction  since  the  local  Mach  number,  Mt.  ,  seen  by  the 
"inner  bourklaiy  layer"  is  certainly  less  than  the  upstream  )&ch 
number,  M,  .  Even  neglecting  the  "inner  boundary  layer", 

should  increase  downstream  of  a  shockvrave  because  the  Mach  number 
external  to  the  boundary  layer,  /Ve  ,  is  less  than  M,  .  However, 
the  existence  of  an  "iruier  boundary  layer"  creates  an  additional 
Increase  in  (since  hence  an  increase  in 

as  seen  from  the  equation 


A  serious  doubt  is  cast  on  the  "inner  boundary  layer”  concept  since 
all  the  data  investigated,  including  the  results  of  the  present  tests, 
show  no  indication  of  an  increase  in  downstream  of  the  peak 
Stanton  number  ratio.  Figures  8U  and  85  are  typical  of  the  local 
recovery  factor  distributions  found  in  this  study.  It  is  seen  that 
downstream  of  station  12.0  for  the  case  of  an  impinging  shock 

is  consistently  lower  than  the  zero  pressure  gradient  .  However, 
in  Figure  85  a  quite  pronounced  peak  in  %  is  found  near  station  9.5* 
This  rise  in  ,  which  starts  at  the  beginning  of  the  pressure 
interaction  a^  then  ends  a  short  distance  later,  may  be  due  to  the 
creation  of  an  "inner  boundary  layer"  which  is  then  destroyed  before 
the  end  of  the  interaction  region  is  reached.  It  is  to  be  noted  that 
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on  page  22  this  *bllp"  In  ^e.  Is  compared  with  similar  results 
obtained  bjr  Brlnlch  employing  two-dimensional  roughness  strips. 

The  results  of  certain  low  speed  experiments  suggest  an  alternative 
to  the  "inner  boundary  layer”  concept  which  gives  a  qualitative 
explanation  to  the  observed  Stanton  number  ratio  rise,  downstream 
of  a  shockwave.  In  the  present  case  the  use  of  subsonic  experiments 
to  gain  insight  into  the  Interaction  of  a  supersonic  flow  field 
with  a  turbulent  boundary  layer  appears  valid.  For  example,  there 
Is  definite  similarity  with  respect  to  local  Stanton  nund>er  ratio, 
s  distribution  between  the  low  speed  flow  disturbed  by 

a  rod  and  the  N.  ■  2.UU  flow  disturbed  by  a  slab  (see  Flgura  82). 

A  fairly  obvious  approach  that  may  be  used  in  an  investigation  of  the 
properties  of  a  reattaching  separated  flow  Is  to  consider  the  flow 
as  being  equivalent  to  a  wall  jet.  Although  the  analogy  may  not  be 
exact.  It  Is  difficult  not  to  believe  that  the  two  phenomena  are 
basically  comparable.  It  Is  pointed  out  by  Meyers  (Reference  19) 
that  the  tests  of  Reference  (21)  show  that  low  speed  wall  jet 
turbulence  levels**  at  the  point  of  maximum  Jet  velocity  are  almost 
four  times  the  maximum  turbulence  level  found  in  a  zero  pressure  gradient 
boundary  layer.  Also  in  Reference  (19)>  Meyers  suggests  that  the 
ratio  of  eddy  diffusivity  for  heat,  ,  and  eddy  diffUslvlty  for 
Biomentum,  *  varies  in  the  manner  ^own  in  the  sketch  below. 


Variation  of  Eddy  Diffusivity  Ratio  (Taken  From  Reference  19) 


*Tbe  term  "turbulence  level”  is  defined  here  to  be  equal  to  ,  where 
■uf  Is  the  root-mean-square  of  the  /y  component  of  the  turbulent  velocity 
fluctuations  and  £/«  is  a  reference  velocity  outside  the  boundary  layer. 
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Therefore,  not  only  a  high  level  of  tnrbulence  le  found  vithin  a 
vail  Jet,  but  apparently  the  turbulent  mechanism  has  been  modified 
eo  that  thei«  is  a  large  increase  in  the  ability  to  transport  heat 
as  eonqjared  with  the  transport  of  momentum.  If  the  wall  Jet  is  truly 
equivalent  to  a  reattaching  separated  flow,  then  it  may  be  seen 
tron  the  sketch  that  high  heat  transfer  rates,  and  yet  low  values  of 
wall  shearing  stress,  are  obtained  downstream  of  ths  reattachment  point. 
In  laminar  sublayer  at  the  wall  where  ths  turbulence  disappears, 
it  can  be  shown  that  the  eddy  diffuslvity  ratio,  A  ,  is  equal  to 
the  raciprocal  of  the  Prandtl  number.  Since  the  recovery  factor 
was  found  to  decrease  downstream  of  the  ^ckwave  for  almost  all 
supersonic  experiments,  it  is  iiq>lied  that  the  Prandtl  number  also 
decreased.  Thus,  if  in  the  foregoing  sketch  the  appropriate  line 
for  flow  downstream  of  a  shockwave  were  drawn  it  would  probably 
intersect  the  A  axis  at  a  slightly  higher  value  than  the  flat 
plate  line. 


C.  MOMENTUM  THICKNESS  AND  HEAT  TRANSFER  COMPARISONS  FOR  NO 
INTERACTION 

The  theory  of  Reshotko  and  Tucker  (Reference  7)  was  employed 
to  calculate  a  theoretical  momentum  thickness,  &  ,  and  heat  transfer 
coefficient,  hr  $  comparison  with  the  present  experimental  data 
in  the  eases  where  no  interaction  was  present. 

The  theory  of  Reference  (7)  applies  to  a  body  with  a  smooth  leading 
edge.  Also,  the  transition  point  in  the  flow  must  be  known  for 
accurate  application  of  the  theory.  In  the  present  calculations 
turbulent  flow  was  assumed  to  begin  at  the  leading  edge  and  the 
surface  was  assumed  to  be  isothermal  from  the  leading  edge  (the 
heat  transfer  data  presented  herein  are  corrected  to  an  isc^ermal 
surface).  Leading  edge  blowii^  was  used  in  some  of  the  M,  ■  2.98 
tests  and  blowing  plus  grit  were  used  in  the  t(.77  tests.  Also, 

the  origin  of  turbulent  flow  was  not  known  thou^  it  is  certain  that 
the  origin  was  not  at  ths  leading  edge.  Therefore,  conditions  were 
present  in  some  of  the  tests  which  cannot  be  theoretically  accounted 
for. 

The  most  important  aspects  of  the  theory  that  can  be  checked  from  the 
present  data  are  the  theoretical  prediction  of  the  magnitudes  of 
hr  and  6  and  the  oi-  dependence  of  hr  •  In  the  present  case  the 
instrumentation  did  not  allow  a  detailed  survey  of  the  boundaiy  layer. 
Therefore,  the  measurements  of  boundary  layer  thickness,  £  ,  and  0 
are  considered  to  be  quite  approximate.  However,  they  are  sufficiently 
accurate  to  give  results  in  the  pressure  interaction  length  correlations 
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(Figures  71  and  72)  idileh  are  consistent  vlth  those  of  other  tests 
wherein  the  boundaiy  layer  vas  earefnily  probed. 

Zero  pressure  gradient  values  of  Bti  were  calculated  from  the  theory 
of  Reference  (7)  for  the  present  fit  «  2.98  and  fi,  •  U.77  data  and 
for  the  tests  from  Reference  (1)  which  are  used  in  this  report. 

The  equations  from  Reference  (7)  for  flat  plate  flow  over  a  smooth 
isothermal  flat  plate,  turbulent  from  the  leading  edge  are  as  follows. 


O.oxs^i 


^r.J 


1 


(¥) 


(1) 


where 

Z. 

'X  ' 

The  values  Hi^^p 


O.IZJ 


'  T^C/to.i  Mt) 


-(2) 


■  "Incompressible"  form  factor  for  a  flat  plate 
*  Distance  from  leading  edge. 

■  1.28,  rr  ■  0.89  and  appropriate  values  for 


7^  (see  Table  III)  were  used  in  the  calculations  of  A,/*}  .  In 

Table  III  the  values  of  A/e  and  the  resulting  from  the 

experimental  measurements  are  compared  with  the  theoretical  values 
of  .  Also  appearing  in  Table  III  are  the  values  of  for  the 
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prasent  tests  calculated  from  the  theory  using  both  the  values  of 
from  the  rake  measurements  and  the  theoretical  values.  The 
theoretical  hri*)  using  the  measured  and  theoretical  Bp.  *8  and 
employing  the  theoretical  d::  dependence  of  'jc  for  both  the 
■saeured  and  theoretical  6^*8  are  shown  on  all  of  the  heat  transfer 
coefficient  figures  for  the  purpose  of  comparing  the  theoretical 
magnitudes  and  ^  dependence  of  hj-  with  the  experimental  results. 

Table  III  shows  that  for  the  present  tests  the  theoretical  ^ 
overestimates  the  experimental  for  M,  ■  2.98  and  1.  $ 

aiv)  underestimates  It  for  H|  «  2.98  and  *  ^.9.  This  ie  to 

be  expected  since  the  delay  In  transition  would  cause  the  measured  Oti 
to  be  smaller  than  the  theoretical  for  the  case  of  no  blowing.  With 
bonndaiy  layer  hloving  more  than  doubled  and  transition 

occurred  further  upstream  than  In  the  case  of  no  blowing.  Because  the 
measured  for  I  ^fe./p  “5.9  lie  above  and  below 

the  respective  theoretical  ,  the  converse  Is  found  for  the  hr 

values  calculated  from  Equation  (2).  Since  varies  as  the  0.268 
power  In  the  h-r  formula,  the  differences  in  hr  using  the 
different  values  of  are  small  and  all  heat  transfer  calculations 
agree  well  with  the  experimental  zero  pressure  gradient  data 


for  M, 


2.98. 


For  M,  ■  li.77  the  grit  and  leading  edge  blowing  which  were  used 
In  every  turbulent  run  caused  the  measured  to  be  over  twice  as 
high  as  the  theoretical  value.  The  increase  in  blowing  from  U.0 
to  9.1  had  little  effect  on  changing  S  and  Q  at  this  Mach 
number.  It  Is  surprising  that  the  theoretical  value  of  gave 
results  for  hr  which  produced  better  agreement  with  the  experimental 
values  than  did  the  measured  value  of  used  In  the  formula. 

Figures  22  through  25  and  1*7  through  1*9  compare  the  theoretical  hr^e 
(using  theoretical  ^  and  measured  with  the  experimental 

values.  It  Is  seen  that  the  theory  of  Reference  (7)  predicts  the 
magnitudes  of  hr  dependence  of  fairly  well  even 

if  is  off  by  more  than  a  factor  of  2.  The  heat  transfer 
coefficient  appears  to  be  insensitive  to  leading  edge  effects,  sudi  > 
as  grit,  unlike  the  pressure  interaction  length  parameters. 

The  best  agreement  between  theoretical  and  measured  found  in 
Table  III  was  obtained  from  the  data  of  Reference  (1).  From  Mach 
nniiA)er  and  Re;ynold8  number  considerations  It  Is  probable  that 
transition  occurred  closer  to  the  leading  edge  in  the  Reference  (1) 
tests  than  in  the  present  tests,  and  since  the  boundary  layer  probe 
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■easarements  were  thorough,  fairly  accurate  values  for  should 

have  been  obtained.  Thus,  the  theory  of  Reference  (7)  is  partially 
verified  by  the  agreement  with  the  measured  0^'s  of  Reference  (1). 

It  aiq>ears  from  the  foregoing  comparisons  that  the  insensitivity  of 
the  heat  transfer  coefficient  to  the  boundary  layer  parameter  06. 
allows  an  estimation  of  ^  to  be  made  fairly  accurately  by  a  theory 
such  as  that  of  Reference  (7)  even  though  the  boundaiy  layer 
parameters  may  be  poorly  estimated. 

&  the  preparation  of  this  report,  whenever  S  or  0  in  a  zero 
pressure  gradient  field  were  required  at  an  ^  location  upstream  or 
downstream  of  where  the  boundary  layer  was  probed,  they  were  adjusted 
from  their  measured  values  by  using  the  theoretical  ^  dependence 
of  from  Reference  (7).  In  Table  III  the  values  of  S  and  ^ 

obtained  in  the  present  tests  from  rake  measurements  at  model  stations 
9.0  and  19.0  are  shown  for  four  cases  in  idiich  the  generator  plate 
was  removed.  As  mentioned  previously  the  values  of  S  and  Q  obtained 
from  the  fixed  rakes  are  not  as  accurate  as  those  obtained  from 
movable  probes  but  they  may  be  used  to  give  an  approximate  value  of 
the  ^  dependence  obtained  in  these  tests.  It  may  be  seen  in  Table  III 
that  the  theoretical  dependence  for  S  anl  6  is  fairly  close 

to  the  average  value  of  found  from  the  Mi  ■  2.98  data.  However, 

the  measurements  at  M,*^  *<77  show  an  almost  negligible  variation  of  S 
and  ^  with  distance.  This  result  is  not  understood.  The  rake 
measurements  of  local  Mach  number  near  the  edge  of  the  boundary  layer 
indicated  an  expansion  large  enough  to  show  up  on  the  pressure 
distributions;  however,  the  recorded  pressures  do  not  indicate  such 
an  expansion. 

Many  calculations  based  on  6,  and  Q,  ,  which  were  determined  from 
the  measured  values  by  using  the  theoretical  dependence, 

were  completed  before  the  negligible  dependence  of  these 
quantities  at  H,"  U. 77 was  noticed.  Since  the  distances  that  the 
measured  S  and  0  were  moved  to  obtain  S’,  and  0,  were  always 
small  (less  than  9^  of  the  value  of  ^  at  which  the  measurements 
were  made  in  all  cases)  the  c  Lculations  were  not  changed. 


D.  CORREUTION  OF  TRANSFORMED  MOMENTUM  THICKNESS  ACROSS  THE 
INTERACTION  REGION 

In  order  to  calculate  the  heat  transfer  coefficient  from  theories 
such  as  that  of  Reference  (7),  knowledge  of  the  variation  of  6 
with  X  is  required.  This  is  because  ©  is  the  only  varying  boundary 
layer  parameter  Involved  in  the  simplified  Reynold's  analogy  and  skin 
friction  formula.  In  the  region  where  the  shock  interacts  with  the 
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botindaiy  layerj  however^  large  gradients  in  the  thermal  and  fltild 
boundar7  layer  parameters  occur  even  for  weak  shocks  as  evidenced 
by  the  present  data  and  the  data  of  References  (1)  and  (2).  Thus, 
the  gradient  lerms  in  the  modified  Reynold's  analogy  of  Cohen, 

Reference  (lU),  become  sufficiently  important  such  that  the  heat 
transfer  depends  on  more  than  O  alone.  Zhtuitively,  it  would  seem 
that  at  some  distance  downstream  of  this  high  gradient  region  the 
heat  transfer  coefficient  would  again  be  given  in  terms  of  the  one 
parameter,  0  .  The  distazice  downstream  of  the  shock  at  which  the 
simplified  Reynold's  analogy  and  skin  friction  formula  can  again 
be  used  to  give  the  heat  transfer  coefficient  is  not  known.  One 
reason  for  the  uncertainty  is  that  downstream  of  the  interaction  >- 
the  velocity  profile  is  filling  out  from  the  distorted  form,  which 
was  caused  by  the  interaction  with  the  ^ockwave,  so  that  Q  is 
changed  at  an  unknown  rate  with  'X  .  Also,  the  present  recovery 
factor,  ^  ,  data  indicate  that  ^  and  the  Reynolds  analogy  factor, 

£  ,  decrease  (assuming  and  for  this  argument) 

to  a  value  downstream  of  the  shock  which  is  lower  than  expected 
from  the  change  in  Pr  due  to  the  change  in  static  temperature. 

This  indicates  a  change  in  the  transport  properties  of  the  boundary 
layer  caused  by  the  interaction.  The  effects  of  this  change  would 
also  (intuitively)  be  expected  to  disappear  somewhere  downstream. 

Effects  such  as  these  prohibit  the  calculation  of  the  heat  transfer 
coefficient  downstream  of  the  interaction  using  the  sin^^le  Reynold's 
analogy  and  skin  friction  formulae  with  downstream  conditions 
determined  from  the  shock  relations. 

Due  to  the  lack  of  experimental  data  and  theoretical  guidelines, 
it  is  difficult  to  arrive  at  any  definite  statements  concerning  the 
downstream  heat  transfer  decay.  However,  an  attempt  was  made  employing 
an  indirect  approach.  Under  the  assumption  that  the  theory  of 
Reference  (7)  applies  at  and  downstream  of  a  "reattachment "  point 
(defined  below),  approximate  values  of  ^  and  .5  at  the  reattachment 
point  may  be  determined  by  utilizing  the  present  heat  transfer  data 
and  the  boundary  layer  profile  data  from  References  (1)  and  (2). 

These  values  may  serve  as  initial  values  for  calculations  of  heat 
transfer  downstream  of  the  reattachraent  point  if  the  0^  dependence 
of  and  S'  can  be  estimated  in  the  theoretical  formula  for 
Stanton  number. 

The  correlation  of  the  transformed  momentum  thickness  ratio, 
across  the  interaction  will  now  be  described. 

The  present  data  show  that  the  peak  heat  transfer  coefficient  occurred 
slightly  upstream  of  the  end  of  the  pressure  interaction.  Therefore, 
even  in  cases  where  a  separation  bubble  existed,  the  peak  beating 
occurred  well  into  the  reattachment  region  so  that  the  ordinaiy 
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boundary  layer  equations  can  be  applied.  For  both  attached  and 
separated  flows  It  Is  assumed  that  the  theory  of  Reference  (7) 
applies  at  the  point  Vhere  the  peak  in  the  hr  distribution  occurs. 
This  assanq>tlan  Implies  that  any  effects  on  hr  due  to  significant 
gradients  in  the  flow  field  existing  at  the  peak  hr  location 
will  be  combined  into  the  values  of  and  obtained  by 

means  of  the  method  described  below.  It  is  implicitly  assumed, 
therefore,  that  the  skin  friction  formula  developed  in  Reference  (7) 
applies  near  reattachment. 

The  formula  for  Stanton  number  from  Reference  (7)  is 


St  ^ 


h  o.u3  [-LSi/ 


s-Ce*) 


.ftztt 


(3) 


where 


St 

% 


“  Incompressible"  form  factor  for  flat  plats  ■  constant 

Applying  Equation  (3)  at  the  beginning  of  the  pressure  interaction, 

/i:,  ,  and  at  where  the  peak  heat  transfer  coefficient,  t 
occurs  gives  the  following  equation  for  i  ratio 

of  the  transformed  momentum  thickness  across  the  heat  transfer 
interaction  region. 


^U.1 


VMi/ 
/  / 


hP 
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Using  reflected  shock  relationships,  the  external  conditions  art  the 
peak  heat  transfer  location  were  detennined  from  plots  of  fU.  and 

versus  pressure  ratio  with 

the  upstream  Kach  number  M,  as  parameter  (see  Figures  75#  76# 

86  and  87).  The  wall  static  pressure  ratio  measured  at  the  peak 
heat  transfer  location  was  used  in  finding  these  external  conditions 
in  each  case.  All  values  used  in  the  calculations  appear  in  Table  H. 
With  h,  and  /r^  known  from  the  data.  Equation  (h)  may  be  used  to 

a.7J 

calculate  the  product,  J  •  Since  Sx  1®  not 

known,  use  was  made  of  shock  wave  boundary  layer  interaction  data 
from  other  sources.  That  is,  boundary  layers  were  probed  through 
tbs  interaction  region  in  References  (1)  and  (2);  the  first  using 
a  flat  plate  and  covering  a  range  of  Mach  numbers,  ,  and 

flow  deflection  angles j  the  second  using  the  tunnel  wall  boundary 
layer  and  one  Mach  number,  one  ,  and  various  shockwave  angles. 

Seven  values  of  were  obtained  from  the  data  of  Reference  (1) 

and  one  value  from  Reference  (2),  The  position  in  the  interaction 
region  at  which  Otkx  was  taken  in  these  two  experiments  was  chosen 
to  correspond  as  closely  as  possible  to  where  the  peak  heat  transfer 
coefficient  would  be  expected  to  occur  in  the  present  tests.  This 
position  was,  for  >  6°,  approximately  255C  of  measured 

upstream  from  .  Since  the  probes  were  not  always  positioned 
idiere  desired,  an  additional  approximation  to  Oicz  was  introduced. 

Uy  comparing  the  ratios  obtained  from  the  direct  measurements 

with  the  values  of  [5*''^']  *  ^  indirectly  from  the 

present  heat  transfer  data,  it  was  noticed  that  the  two  sets  of 
values  could  be  brought  into  general  agreement  If  the  value  of 
rs,/  13.73 

[_  was  taken  as  U  in  all  cases.  This  means  that 

■  0.690  S',  ,  where  -S',  ~  •  0.791.  This  decrease 

in  S  from  S,  to  S'^.  is  in  line  with  the  observed  decrease  in 
the  experimental  recovery  factors  from  7?/  to  .No  comparison 
with  directly  measured  values  was  available  for  the 
ratios  calculated  from  the  present  M,  *  U.77  data. 
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The  mlues  of  are  listed  in  Table  II  and  are  plotted 

as  functions  of  H,  vith  tl»  parameters  (C  and  in  Figure  88. 

It  is  seen  in  Figure  88  that  the  values  fall  into  three 


rather  distinct  bands  even  though  the  error  involved  in  detenaining 
the  ratio  from  the  heat  transfer  data  is  quite  large  (due  to  the  fact 

I  I  and  the  error  in  ht  and  are  by  nature 


large).  Table  II  also  indicates  whether  the  flow  is  attached. 
Incipient,  or  separated  for  each  of  the  tests.  It  nay  be  seen  from 
Figure  88  that  Ibch  number  and  boundary  layer  blowing  have  a  small 

effect  on  but  there  is  no  discemable  trend  due  to 


The  condition  of  the  flow,  whether  separated,  incipient  or  attached, 
was  determined  by  applying  Kuehn's  criterion  (Reference  6)  idiich  was 
discussed  in  Section  III  A  1.  .  Kuehn's  results  were  used  to  define 
the  condition  of  the  flow  for  the  data  of  References  (1),  (2)  and  (6) 
and  the  present  H,  *  2.98  data.  Since  Kuehn's  data  did  not  extend 
to  M,  •  I4.77,  this  case  had  to  be  Judged  from  the  pressure  data 
taken.  In  the  present  experiment  the  spacing  of  the  pressure  taps 
allowed  the  Incipient  separation  conditlm  to  be  determined  only  to 
within  about  •  t  !<>.  It  is  concluded  that  the  oC  ■  9.7°  and 
t>C  ■  11°  tests  at  M,  ■  It. 77  correspond  to  Just  below  incipient 
and  Just  above  incipient  separation  respectively. 
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V.  DBSICaf  PBOCEDDRE 


In  this  section  a  step-ly-step  procedure  is  outlined  for  calculating 
the  isothermal  heat  transfer  coefficient  distribution  idiLch  result 
from  an  tnqpinging  oblique  shockwave  interacting  with  a  turbulent 
boundary  layer  on  a  flat  plate.  The  flow  field  is  considered  to  be 
two-dimensional  in  the  planar  sense.  The  essential  features  of  the 
calculation  are  the  estimation  of:  (1)  the  length  of  the  pressure 
interaction  region,  4^  ,  as  determined  by  the  static  pressure 
distribution  on  the  flat  plate  and  (2)  the  isothermal  heat  transfer 
coefficient,  ^  ,  distribution  from  tbs  beginning  of  the  pressure 
interaction  to  a  point  a  short  distance  downstream  of  the  end  of  the 
Interaction.  Correlations  of  the  interaction  length  were  discussed 
in  Section  I7-A  idiich  Indicated  that  the  method  of  inducing  transition 
had  a  strong  influence  on  the  correlation.  However,  in  this  section 
A/)e  is  determined  from  correlations  which  are  more  appropriate  for 
use  in  a  design  calculation.  The  calculation  of  the  teat  transfer 
coefficient  distribution  resulting  from  the  shockwave  turbulent 
boundary  layer  Interaction  was  obtained  primarily  from  the  tentative 
relationships  discussed  in  Section  IV-B. 

The  calculation  procedure  requires  that  the  following  quantities 
are  known: 

1.  The  Kacb  number  upstream  of  the  interaction,  M| 

2.  The  boundary  layer  thlckxiess  Immediately  upstream  of 
the  interaction,  S,  ,  hn  equation  for  estimating  S, 
can  be  obtained  from  Reference  (1): 


where  U,  and  Vi  are  the  velocity  and  the  kinematic  viscosity, 
respectively,  in  the  Inviscid  flow  upstream  of  the  interaction. 

L.  is  the  distance  from  the  point  of  bouxidary  layer  transition  to 
the  beginning  of  the  interaction,  and  A/  is  the  reciprocal  of  the 
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expomnt  in  tbe  boundary  layer  power  profile  lav, 


/ 


3*  The  point  on  the  flat  plate,  ^xk  »  which  corresponds  to 
the  intersection  of  the  oblique  shockwave  with  the  edge  of 
the  boundary  layer  if  there  were  no  Interaction. 

This  is  illustrated  in  Figure  92. 

Ii.  The  Reynolds  number  per  unit  length  upstream  of  the 
interaction,  . 

5.  The  static  pressure  ratio  across  the  oblique  reflected 

shockwave  in  an  inviscid  flow  field,  (see  Figure  89). 

6.  The  isothermal  heat  transfer  coefficient  upstream  of  the 
interaction,  h,  •  A  theoretical  prediction  of  f), 

by  the  theory  of  Reference  (7)  can  be  obtained  by 
enploying  Equations  (1)  and  (2)  (which  are  valid  for  an 
isothermal  wall)  in  Section  I7-C. 

The  calculation  procedure  given  below  pertains  only  to  a  wall  tdiich 
is  isothermal  far  upstream  of  the  interaction. 


Step-by-Step  Calculation  Procedure 

Step  1:  Calculate  the  pressure  interaction  length, 

from  either  Figure  90  or  91.  Figure  90  will  give  greater 
accuracy  if  it  is  known  whether  the  boundary  layer 
transition  is  natural  or  forced. 


Step  2;  Locate  the  interaction  length  in  relation  to  by 

calculating  from  Figure  92. 

Step  3!  From  the  known  valu?  of  e<  and  find  the  mass 
flow  ratio,  across  the  reflected  oblique 

shock  from  Figure  93  or  83. 


NORTH  AMERICAN  AVIATION 

COIMMtUS  »IVI$IOM 
COtUMtVS  l«.  OHIO 


INC. 

HA  62H-79^ 


Step  U: 


The  peak  value  of  the  Stanton  nunber  ratio, 
can  also  be  determined  with  the  aid  of  Figure  93« 

If  the  combination  of  and  gives  a  point  on 
the  curve  for  a  given  M,  to  the  left  of  the  dashed 
enrre  in  Figure  93»  then  «  1,6  or  1.8 


(Hote:  K|  ■  1.8  is  the  highest  value  of  found 

experimentally  for  all  ^  average 


peak  value  of  all  the  present  data,  see  Sketch  on  page  31). 
If  0^  and  t1,  locate  a  point  to  the  right  of  the  dashed 
curve  in  Figure  93,  then 


-1%)^  J 


Step  5: 


For  a  given  Kach  number.  Mi 


the  pressure  ratio 


corresponds  to  the  value  of  ^  at  the  intersection 

of  the  dashed  line  with  a  solid  line.^F^^  j  can  be 


found  directly  from  Figure  83  or  by  determining  the 
pressure  ratio  in  Figure  89  which  corresponds  to  0^15 
and  M,  . 

Calculate  from: 


I 


ery  1 


Step  6:  Sketch  the  /if  distribution  from  oc,  to  the  peak  value  /i^ 
at  according  to  Figure  9U  as  follows:  Extend  the 
upstream  value  of  h  15^  of  downstream  of  Ac,  j 

locate  'Xi.  25^  of  upstream  of  ;  connect  the 

two  points,  and  »  by  a  straight  line. 

Then  extend  the  peak  value,  ht  »  downstream  to  /Xj  where 
the  downstream  decay  will  begin. 
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step  7:  Rram  Figure  9$  determine  the  downstream  decay  of  the  heat 
transfer  coefficient  beginning  at  .  The  rate  of  decay 
was  determined  from  Figure  8l.  The  isothermal  beat 
transfer  coefficient  downstream  of  X5  is  calculated 
from  the  formula, 

1.^=  hjlh] _ ^ - 

/  Ft**  Secx  •R 

Append  the  downstream  decay  of  Vi  ^  to  the  K.|.  distribution 
obtained  in  Step  6  beginning  at  x ^  . 

Step  8:  Calculate  the  recovery  tenqierature 

with  ^  ■0.89  where  Tt  represents  the  freestream 
stagnation  temperature.  Use  this  ell  the  way 
through  and  downstream  of  the  interaction. 

Step  9:  Calculate 

using  the  k.|.  distribution  obtained. 


A.  EXAMPLE  CALCUUTION 

fiiown  conditions:  M|  ■  U.77»  ■  0.3li",  ■  3.32  x  lO^/in,, 

10.7"  from  leading  edge,  ■  7.25  (foraC*  11°) 

Step  1:  From  Figure  89  or  76,  of  ■  11®  and  M,  ■  U.77,  it  is  inferred  that 

■  8.25.  However,  due  to  the  expansion  from  the  shock 
generator  trailing  edge  the  peak  pressure  ratio  is  reduced  to 
a  value  of  7.25  (see  Figure  U3).  Use  Figure  91  (assuming  type 
of  transition,  forced  or  natural,  is  not  known)  to  obtain  AX  . 


-  ■  1.22  Figure  91  gives,  ■  lU.U 


So,  AX  ■  U.9  in. 
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Step  2:  With 


V.-' 


■  I1I46  Figure  92  givee. 


.XtK-X|  _0  33,^y 

So,  >^|  ■  10.7  -  (0.325  at  U.9)  ■  9.1  in,  from  leading  edge, 
X^r  X,  +  AX  -  Ui.O  in. 


Steps  3  and  It:  Figures  93  and  89  or  76  gives  the  peak  value 

f  1.8  conservative  experimental  value  (jPU)^ 

""  ”1^1.6  average  eacperlmental  value  **** 

Step  The  theoretical  value  of  It,  (from  Table  III)  for  this  test 
run  at  X  •  9.0"  vas  h,  ■  2,2  x  10"3  (BTU/ft^  x  sec  x  °R). 

Therefore,  “  2.2  x  10“3  x  ^  x  3.5 

-  13.9  X  10"^  Is  s  conservative  estimate 

of  the  peak  heat  transfer  coefficient  and  hj.  “  2.2  x  10~3 

X  1,6  X  3.5  ■  12.3  X  10"^  (■—s — — r; — )  represents  an 
-  ft‘  X  sec  X  ®R 

average  value  which  shows  better  agreement  with  the  present 
dfata. 

Step  6:  See  Figure  96  for  the  sketch  of  the  Hi*  distribution  and 
caqurison  with  the  experimental  distribution. 

Step  7:  In  Figure  96  Kx  downstream  of  Xj  was  calculated  using 
Figure  95  From  the  formula, 


l-T=  >< 


The  coordinate  X  was 


obtained  from  the  dimensionless  coordinate  in  Figure  95  by 
using  the  from  Step  2  and  X)  from  Step  3.  Note  that 
the  impingement  of  the  expansion  wave  from  the  trailing  edge 
of  the  shock  generator  plate  tends  to  reduce  K<|-  at  a 
faster  rate  than  would  occur  if  no  expansion  were  present 
so  that  comparison  of  the  experimental  to  the  estimated 
is  not  possible  beyond  the  expansion  wave  impingement  point. 
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?I.  CONCLPDIMO  REKURKS 


A.  REVIEW  OP  skruficaiit  results 

The  mast  elgnif leant  results  found  in  this  study  are: 

1.  The  isothermal  heat  transfer  coefficient,  hj.  , 
on  the  undisturbed  test  plate  nay  be 

accurately  predicted  by  the  theory  of  Reference  (7), 
eren  though  the  use  of  grit  and  air  Injection  into 
the  boundary  layer  gives  measured  values  of 
momentum  thickness,  &  ,  considerably  different 
fromd  predicted  by  the  theory  of  Reference  (7)* 

The  insensitivity  ot  to  0  predicted  by 
theory  appears  valid. 

2.  Ih  Figures  71  and  72  plots  are  given  of  pressure 
interaction  length,  ,  as  a  function  of  other 
parameters  which  characterize  the  flow.  The  data 
correlate  along  straight  lines  with  each  line  Indicative 
of  the  method  of  inducing  transition.  From  these  figures 
the  following  observations  are  noted; 

(a)  Artificial  methods  of  inducing  transition  yield 
a  boundary  layer  with  a  turbulent  structure 
ihich  is  different  from  the  mechanism  in 
"naturally"  turbulent  flow. 

'  (b)  Moderate  amounts  of  air  injection  into  a 

boundary  Iryer  do  not  alter  its  turbulent 
structure.  Thus,  if  air  injection  or  blowing 
thickens  a  turbulent  boundary  layer,  then  its 
subsequent  behavior  will  be  as  if  the  thickening 
was  the  result  of  a  longer  run  of  an  undisturbed, 
zero  pressure  gradient,  turbulent  flow. 

3*  The  adiabatic  wall  tesq^erature,  ,  and  the  effective 
recoveiy  temperature.  Tit  »  approximately,  but  pot 
exactly,  the  same  quantity.  (Note:7(  is  the 
temperature  at  which  the  local  heat  transfer  rate  is  zero. 
The  average  value  of  TH  and  in  a  shockwave  boundary 
l«yer  interaction  region  are  nearly  equal  but  the 
variations  in  Tl  are  between  wider  limits  than  is  found 
for  . 
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li«  At  Mach  5  for  the  eases  of  shock  impingement  on  the 
test  plate^  a  "blip*  vas  found  in  the  pecoveiy 
factor,  ,  distribution.  (Note:  X  is  proportional 
to  ).  This  blip  is  similar  to  the  results 
obtained  in  Reference  (11)  with  two-dimensional 
roughness  elements. 

5.  The  peak  in  the  heat  transfer  coefficient  occurs 
iq}stream  of  the  end  of  the  pressure  interaction 
region.  The  sero  pressure  gradient  value  of  heat 
transfer  coefficient  is  maintained  a  short  distance 
downstream  of  the  beginning  of  the  pressure  inter¬ 
action.  The  heat  transfer  coefficient  then  rises 
monotonically  (except  perhaps  over  a  very  small 
distance)  to  the  peak  value. 


6.  The  Stanton  number  ratio,  ,  is  defined  as  the 

ratio  of  the  Stanton  number  based  on  local  flow 
conditions  external  to  the  boundary  layer  to  the 
zero  pressure  gradient  Stanton  number.  The  local 
conditions  external  to  the  boundary  layer  may  be 
estimated  from  the  static  pressure  distribution 
on  the  test  plate.  Plots  have  been  made  of  Sts./ 


versus 


A')' 


,  where  /)C.  is  the  station 


at  trtiich  the  pressure  interaction  begins  a«j 
is  the  length  of  the  pressure  Interaction,  All  of 
the  appropriate  data  obtained  in  the  present  tests 
show  qualitatively  similar  distributions  downstream 
•  immediately  downstream  of 
decreases  to  a  valuer  less  than  unity.  Then 
the  Stanton  number  ratio  Increases  rapidly  and 
usually  reaches  a  peak  in  the  vicinity  of 

-  ■  ‘T?.  For  all  of  the  present  data  the 


peak  in  Vyt,  is  between  1.5  and  1.8.  The  decay 
of  downstream  of  the  peak  is  fairly  slow 

with  remaining  well  above  unity  for  those 

tests  in  which  the  expansion  from  the  trailing  edge 
of  the  shock  generator  plate  did  not  have  a  stroxig 
influence  on  the  test  plate  static  pressure 
distribution. 
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7*  The  reeulie  x*eported  in  the  literature  of  heat 
transfer  tests  on  reattaching  separated  boundary 
layers  in  a  two-dimensional  flow  field  reveal  an  interesting 
onalitative  similarity  to  the  present  experiments;  vis., 

(1)  in  the  separated  zone  the  Stanton  nukber  ratio  is 
^low  unity  and  then  peaks  in  the  region  of  reattachmsnt, 
and  (2)  the  recovery  factor  based  on  local  conditions, 

,  in  and  downstream  of  the  separated  region  is  less 
than  the  zero  pressure  gradient  value. 

8.  Based  on  the  present  data  and  the  few  published  heat 
transfer  measurements  in  the  Interaction  region  of  a 
shockwave  and  a  turbulent  boundary  layer,  it  is  suggested 
that  the  Stanton  number  ratio  at  the  peak,  has 

the  following  form,  ' 


5t^ 

St, 


K, 


where  is  the  static  pressure  ratio  across  the 
interaction  for  a  given  Mach  number,  is 


the  pressure  ratio  which  corresponds  to  the  "elbow" 
in  the  curve  of  mass  flow  ratio  versus  pressure  ratio 
(e.g.,  see  Figure  83).  Ki  and  Kj.  are  approximately 
constant  with  l.U  to  1.8  and  Kv^  0.1U5. 


9*  A  correlation  of  the  ratio  of  the  transfonned  momentum 
thickness  across  the  heat  transfer  interaction  region, 
,  appears  feasible.  The  values  of  the 

ratio  fall  into  roughly  three  groups  characterized 
by  attached.  Incipient  and  separated  flow  conditions. 
The  indirect  calculation  of  6tva. i^m  the 
present  data  utilizing  many  assumptions  Indicates 
a  decrease  in  the  Reynold's  analogy  factor  across  the 
interaction  region.  This  is  In  qualitative  agreement 
with  the  observed  decrease  in  the  experimental  local 
recoveiy  factor  across  the  interaction  region. 


50 
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B.  RBCOMHElinilTlDNS  FOR  FURTHER  SlUET 

The  dlecueelon  in  the  preceding  sections  indicate  that  certain  aspects 
concerning  the  interaction  of  a  turbulent  boundary  layer  with  a 
shockwave  require  further  study*  First,  in  Section  IV-A  it  was 
pointed  out  that  gross  effects  on  the  pressure  interaction  length 
correlations  are  due  to  the  nelhod  of  inducing  transition  (natural 
or  artificial)  and  that  different  boundary  layer  trlpp^  techniques 
■ay  Induce  diverse  results.  However,  boundary  layer  blowing  seemed 
to  effectively  increase  the  boundary  layer  and  momentum  thicknesses 
without  introducing  spurious  effects  into  the  correlations*  The 
pressures  obtained  with  the  pitot  rakes  showed  that  no  gross 
distortions  of  the  velocity  profile  were  created  by  employing  grit 
and/or  blowing*  This  lnq)lies  that  the  tine-average  pressures  through 
the  boundary  layer  are  not  strongly  affected*  Since  the  boundary 
layer  rake  could  not  measure  the  laminar  sublayer,  the  fraregoing 
observations  might  be  an  Indication  that  the  grit  introduced  a  long 
lasting  disturbance  into  the  sublayer  which  in  turn  exeii^ed  a  large 
influence  on  the  development  of  the  boundary  layer  when  disturbed 
by  the  incident  shock*  Therefore,  it  would  be  desirable  to  conduct  a 
shockwave  turbulent  boundary  layer  test  in  which  various  trips  are 
employed  under  identical  conditions*  Accurate  measurements  should  be 
made  of  the  Interaction  length  and  the  boundary  layer  velocity  profiles, 
especially  the  ^bsonlc  portion,  at  various  stations  beginning 
immediately  downstream  of  the  trip* 

In  the  present  tests  a  qualitative  difference  in  the  local  recovery 
factor  distributions  be'^een  Hach  nunbers  3  and  5  was  found*  At 
Kach  3  the  recovery  factor,  decreased  somewhat  through  and  downstream 
of  the  interaction  region,  but  at  H«  •  $  the  local  recovery  factor 
increased  to  a  high  value  near  the  beginning  of  the  interaction  region 
and  then  decreased  to  a  value  dowi:stream  of  the  interaction  which  was 
lower  than  the  upstream  value*  The  reasons  for  this  rather  profound 
difference  between  the  H,  ~  3  and  H,  ■  5  data  are  not  understood* 
Additional  experiirents  should  be  ccmducted  in  which  adiabatic  wall 
teiig>eratures  are  measured*  Perhaps  tenperature  measurements  made  in 
oonjunction  with  the  tests  suggested  in  the  preceding  paragraph  would 
isolate  the  factors  which  are  the  primary  cause  of  the  effect* 

The  i>resent  data  showed  that  the  Stanton  number  is  high  at  aid  down¬ 
stream  of  the  point  where  a  separated  turbulent  boundary  layer  reattaches* 
However,  Qadd  in  Reference  (8)  and  Seban  in  Reference  (18)  indicate 
that  the  skin  friction  is  low  in  this  region*  It  is  probable  that 
low  values  of  skin  friction  would  also  be  found  downstream  of  a  diock- 
wave  tuztulent  boundary  layer  interaction  even  if  no  boundary  layer 
separation  occurred*  As  discussed  in  Section  IV-B  the  physical  reason 
underlying  the  increase  in  Stanton  number  and  decrease  in  skin  friction 
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eoeffloient  across  the  Interaction  region  nay  be  the  result  of  a 
change  In  the  ^  distribution  of  the  ratio  of  the  heat  to 
■oasntun  eddy  diffostrity,  A( f)  •  It  vas  also  suggested  in 
Section  IV-B  that  accon;>anylng  the  change  in  ia  an 

overall  increase  in  the  level  of  turbulence  in  tne  boundary  layer 
throni^  the  interaction.  It  has  been  shown  experimentally  that  for 
increasing  distance  in  the  downstream  directionf  the  skin  friction 
coefficient  increases  and  the  Stanton  number  decreases  slowly  fron 
their  values  at  the  end  of  the  interaction.  Therefore^  it  is 
eaq>ected  that  both  the  level  of  turbulence  and  vill  eventually 

rstam  to  the  usual  zero  pressure  gradient  values.^  Althou^  the 
measurenents  would  be  tedious  and  eaqpensive,  it  is  felt  that  it  would 
be  worthwhile  to  detenaine  the  variation  in  turbulence  and 
in  the  streamrise  direction  through  and  far  downstream  of  a  shock- 
wave  turbulent  boundary  layar  interaction.  The  physical  insight 
gained  from  such  detailed  measurements  would  probably  allow  a 
theoretical  model  of  the  interaction  to  be  constructed  which  would 
hare  quite  general  application. 

Certain  interesting  properties  of  the  boundary  layer  downstream  of 
the  interaction  region  could  be  determined  from  measurements  which 
are  considerably  siDq)ler  and  less  expensive  to  obtain  than  those 
suggested  in  the  preceding  paragraph.  In  particular  the  measurements 
would  consist  of  accurately  determining  the  heat  transfer  rates  and 
the  boundary  layer  velocity  profile  at  numerous  stations  downstream 
of  the  interaction  region.  From  this  information  it  would  be  possible 
to  determine  tinder  what  conditions,  if  any,  that  the  couparatively 
simple  relationships  between  heat  transfer  coefficient,  momentum 
thickness,  Reynolds  analogy  factor,  and  distance  in  the  streamwise 
direction  given  bf  Reshotko  and  Tucker  in  Reference  (7)  are  valid 
downstream  of  a  shockwave  turbulent  boundary  layer  Interaction. 

Of  course,  it  would  be  usefhl  to  extend  the  range  of  the  present 
types  of  measurements  to  other  I&di  numbers  and  hi^er  pressure 
ratios  along  with  variations  in  Reynolds  number.  Extending  the 
range  of  pressure  ratio  across  the  interaction  region  is  Inportant 
for  the  verification  or  invalidation  of  the  assunption  utilized 
in  the  peak  Stanton  number  formula  in  Section  I7-B  whi(di  considers 
the  peak  Stanton  number  to  be  a  function  of  shock  strength  at  a 
given  Mach  number. 
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m.  HoiATiDir 


q;>eed  of  sotuad 

,  pru 

Q>6elfle  heat  at  constant  pressure,  0.2U  for  air 

aoceleratlon  due  to  gravit7,  32.2  ft/sec^ 


boundary  layer  foxn  factor  (%) 
othervlse  defined 

heat  transfer  coefficient 


exe^  where 


heat  transfer  coefficient  on  a  wall  which  is 
Isothermal  along  the  entire  run  of  the  boundary  layer 

Kach  number 


pressure 

pressure  in  test  plate  plenum  chamber 
Prandtl  number 

heat  transfer  rate  per  unit  area 

Reynolds  number  per  inch  of  free  stream  flow 

Reynolds  number  based  on  S' 

Reynolds  nusRjer  based  on  ^ 

Stanton  nunber 
Roynold's  analogy  factor 
teiqwratnre 

local  static  tenperature  external  to  the  bounda^ 
layer, 

"effective"  recovery  tenperature,  see  Figure  16 
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SUBSCRIPTS 


7*  stagnation  oondltions 

^  eondition  at  edge  of  boundary-  layer 

^  traneforned  to  inempreaaible  plane 

£  related  Inconpresaible  value 

evaluated  at  reference  condition 
o  Oree-etream  value  ig>8treaa  of  test  plate 

P  ^  location  at  which  boundary  layer  was  probed 
refera  to  incident  shock 
AW  adiabatic  wall 

i  eondition  at  beginning  of  pressus  interaction  region 

or  the  condition  for  the  case  of  no  interaction 

a  condition  at  ^  location  where  peak  heat  transfer 
ooefficlent  or  peak  Stanton  number  ratio  occurs 

3  In  referring  to  e3q>erliaBntally  determined  values^ 

tile  subscript  represents  quantities  evaluated  at  the 
location  of  the  experimental  peak  in  the  static  pressure 
distribution  on  the  test  plate.  In  cases  pertaining 
to  a  reflected  shockwave  in  invlscid  flow  tiie  subscripvt 
refers  to  the  tiieoretical  static  conditions  downstzuan 
of  the  reflection. 
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If  It  is  assuiiBd  that  tha  Stanton  nninber,  ,  oorreq;>oi]ding  to 
the  peak  beat  transfer  coefficient,  and  the  location  of  the  peak,  $ 
ean  be  estimated  bgr  a  method  such  as  the  procedure  outlined  in 
Section  T,  then  the  transformed  momentum  thickness,  d»ix  ,  at  4c.x 
can  be  obtained  from  the  correlation  ehovn  in  Figure  88*  According 
to  the  assuiqations  stated  in  SectiLcm  I7>D  the  Stanton  number  Tarlation 
doimstream  of  /Hx  can  be  calculated  from  a  knowledge  of 
as  long  as  the  ^  dependence  of  Qu.  and  the  Re(7nolds  analogy-  factor, 
^  ,  can  be  obtained.  Eqfuation  (3)  in  Section  IV>D  evaluated  at 
4C»  ■1C,  is 


(5) 


where 


r  ^  1  e^  .f/.S6/  Ht,fr  ] 


■7M 


In  general,  conditions  external  to  the  boundary  layer  at  will 
not  be  known  since  the  present  da-ta  shows  that  the  peak  heating 
occurs  somewhat  before  the  wall  pressure  gradient  -vanishes  at  • 

Calculations  eoploylng  Equation  ($)  (using  theoretical  downstream 
e<niditi(»is  for  f  ,  Sx  *  0.690  x  0.791  ■■  0.^6,  and  Q\x 
determined  by  the  correlation  in  Figure  88)  gave  values  for 
wliidh  were  within  20^  of  tte  experimental  values  obtained  in 
the  uesent  tests  for  all  data  runs  except  the  •  6^,  >V/  ■  2.98, 
and  *  1  nn  where  the  calculated  was  $2%  hig^r  than  the 

eaperlmsntal  value.  Although  the  estimate  of  Am.  Equation  ($) 
using  the  theoretical  downstream  conditions  gave  a  poor  estimate 
cmljr  in  the  one  case  mentioned  above,  it  is  felt  that  -the  simpler 
Stan^  nmnber  correlation  of  Section  T7-B  should  yield  more 
consistent  results  In  estimating  I  Equation  (b)  can  iSen  be 
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vsed  in  the  Ibllawlng  way  to  obtain  Stanton  muriber  or  haat  transfer 
eoeffieient  downstream  of  the  peak  heating  point*  Rewriting 
Eiiiiation  (5)  as} 


where  'b*  is  a  constant  eraloated  in  Step  5  below  and  p  is 
evaloated  using  the  theoretical  conditions  in  the  inviscld  flow 
downstream  of  the  interactionjand 

The  proposed  method  is  as  follows: 

1*  Find  Stx  and  ^7  the  procedures  given 

in  Section  V* 

2*  For  the  given  incident  shock  strength  and  M, 
from  the  dt«.  correlation  In  Figure  88* 

3*  Calculate  Oix.t  from  a  theory  such  as  that  of 

!»•  Calculate  ^ 

5*  Evaluate  Equation  (6)  at  using  5L  *  o.tfoi*  ^  £>,yv<6 

BUx  and  Six.  to  obtain  t  ,  't  is 

the  constant  which  brings  both  sides  of  Equation  (6)  into 
agreement  at  . 

6*  With  't*  known^  use  Equation  (6)  with  SC')c}  j 
to  calculate  StC'y)  for  ' 

At  present  for  are  not  known.  One  would 

intuitively  e:q»ct  S(*]  to  change  from  the  value  to  its  value 

for  a  normal  turbulent  boundary  layer  at  some  distance 

downstream  where  the  velocity  profiles  have  filled  out  to  normal 
turbulent  profiles.  Concerning  OxiC^)*  present  boundary  layer 
rake  data  at  Mf  •  2.98  and  the  probe  data  from  Reference  (1)  used 
in  this  report  were  examined  briefly  to  determine  the  order  of 
magnitude  of  the  /y^  d^endence  cf  .  The  distance  in  the 


,  find 

Reference  (?)• 
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strMumrlse  direction  between  and  the  neareat  point  at  which  the 
boondary  layer  was  probed  was  approximately  2"  In  the  present  testa 
and  l.R”  In  Reference  (!)•  The  data  indicated  that  for  an  ecpaatlon 
of  the  fora 

&ti  ‘  ( fj 


J/  was  wllhln  the  limits  Since  the  ^  separation  need  was 

of  the  order  of  A'X  f  V  conU  be  greater  than  3  near  •  It  Is 
eaqacted  that  A/^*)  will  decrease  with  increasing  Af'  from  its 
large  Initial  value  at  as  the  velocity  profile  Is  restored. 

If  A/C'V.)  could  be  determined,  S(ic)  could  probably  be  given  an 
average  value  and  StCif^  for  AO'Ki  could  be  calculated  from 
Equation  (6)  for  cosparison  with  experimental  values.  The  present 
data  give  only  limited  information  concerning  the  downstream  decay  of 
heat  transfer  since  expansion  waves  from  the  trailing  edge  of  the 
generator  plate  usually  Interfered  with  the  heat  transfer  measurements 
a  short  distance  downstream  from  .  In  the  present  tests  the 
greatest  distance  between  and  the  point  where  the  eoqpanslon 
wave  from  the , generator  plate  struck  the  boundary  layer  was  about 
10  upstream  boundary  layer  thicknesses. 
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LET  OF  MTA  HINS  AND  RDM  RdRAMETERS 


Number  M. 

o 


0^ (Degrees) 


1 

2 

3 

U 

$ 

6 

7 

8 

9 

10 

11 

12 

13 

lU 

1$ 

16 

17 

18 

19 

20 

21 

22 

23 

2k 

2$ 


2.95  2.98 


f  t 

5.02  k.77 


V  V 

2.95  2.98 


5.02  U.77 


V 


1.0 

6 

5.9 

6 

1.0 

9.5 

5.9 

9.5 

1.0 

12 

2.0 

12 

5.9 

12 

ii.05 

7.1* 

3.9 

9.7 

8.U 

9.7 

li.05 

11 

6.96 

11 

9.5 

11 

1».9 

12 

9.6 

12 

1.0 

0 

5.9 

0 

1.0 

generator  out 

2.0 

generator  out 

3.7 

0 

1».2 

0 

1(.0 

generator  out 

9.1 

generator  out 

1.5 

10 

l.U 

11 
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Lese  values  of  peak  Stanton  number  ratio  differ  somewhat  from  those  found  in 
^ures  77'80.  This  is  due  to  a  somewhat  arbitraxy  method  of  fairing  the  pressure 
id  heat  transfer  data  to  obtain  the  numbers  used  in  this  table. 
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Distance  from  coordinate  origin  at  idilch  boandarjr  layer  opi 
Interaction  was  probed;  ^  i**  general. 
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whez*e  a  probe  neasorement  of  had  been  made  so  the  approximation  to  de. 
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